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ABSTRACT 

We present echelle long-slit optical spectra of a sample of objects evolving off 
the Asymptotic Giant Branch (AGB), most of them in the pre-planetary neb- 
ula (pPN) phase, obtained with the ESI and MIKE spectrographs at the 10 m 
Keck II and 6.5 m Magellan-I telescopes, respectively. The total wavelength range 
covered with ESI (MIKE) is ~3900 to 10900A (~3600 to 7200A). In this paper, 
we focus our analysis mainly on the Ha profiles. Prominent Ha emission is de- 
tected in half of the objects, most of which show broad Ha wings (with total 
widths of up to ~4000km s _1 ). In the majority of the Ha-emission sources, fast, 
post-AGB winds are revealed by P-Cygni profiles. In ~37% of the objects Ha 
is observed in absorption. In almost all cases, the absorption profile is partially 
filled with emission, leading to complex, structured profiles that are interpreted 
as an indication of incipient post-AGB mass-loss. The rest of the objects (~13%) 
are Ha non-detections. We investigate correlations between the Ha profile and 
different stellar and envelope parameters. All sources in which Ha is seen mainly 
in absorption have F-G type central stars, whereas sources with intense Ha emis- 
sion span a larger range of spectral types from O to G, with a relative maximum 
around B, and also including very late C types. Shocks may be an important 
excitation/ ionization agent of the close stellar surroundings for objects with late 
type central stars. Sources with pure emission or P Cygni Ha profiles have larger 
J — K color excess than objects with Ha mainly in absorption, which suggests 
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the presence of warm dust near the star in the former. The two classes of pro- 
file sources also segregate in the IRAS color-color diagram in a way that intense 
Ha-emitters have dust grains with a larger range of temperatures. Spectral clas- 
sification of the central stars in our sample is presented. For a subsample (13 
objects), the stellar luminosity has been derived from the analysis of the O 1 7771- 
5A infrared triplet. The location in the HR diagram of most of these targets, 
which represent ~30 % of the whole sample, is consistent with relatively high final 
(and, presumably, initial) masses in the range M f ~0.6-0.9M Q (Mi~3-8M ). 

Subject headings: stars: AGB and post-AGB, stars: mass loss, circumstellar 
matter, ISM: jets and outflows, planetary nebulae: general 



1. Introduction 

Intermediate mass stars (~1-8M Q ) evolve from the Asymptotic Giant Branch (AGB) 
to the Planetary Nebula (PN) phase through a short-lived (~ 10 3 yr) and fascinating evolu- 
tionary stage designated as the post-AGB (pAGB) or pre-planetary nebula (pPN) phase. At 
some point in the late- AGB or early pAGB stage, a process (or processes) becomes operative 
that accelerates and imposes severe asymmetries upon the slow, spherical AGB winds: the 
spherical, slowly expanding (V^ xp ~15km s -1 ) AGB circumstellar envelope (CSE) becomes 
a PN with clear departures from sphericity and fast (> 100km s _1 ) outflows directed along 
one or more axis. Although there is no consensus yet for what causes this spectacular meta- 
morphosis, fast jet-like winds ha ye been hypothesized to play an important role (see e.g. the 



review paper on PNs shaping by lBalick fc Frank! |2002| ). These outflows carve out an imprint 



within th e AGB CSE producing and shaping the fast, bipolar lobes observed in most pPNs 



and PNs (jSahai &: Traugerlll998l ). The mechanism that powers and collimates pAGB jets is 



a fundamental issue on stellar evolution that remains a mystery. 

Optical spectroscopic observations of pPNs and PNs have allowed considerable advances 
in our understanding of pAGB evolution and, in particular, are very useful for probing the 
interaction between pAGB winds and the CSE formed in the previous AGB phase. Our 
current (very limited) knowledge of pAGB evolution and, more particularly, of pAGB winds 
is derived mainly in two ways. The first is indirect, that is, based on the effects of the 
pAGB winds on the AGB CSEs. Many pPNs/PNs show extended lobes with, often, bow- 
shaped features at their tips that are visible through optical recombination and forbidden 
emission lines. Spectroscopic observations have been crucial for understanding the origin of 
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these regions, which are found to be excited by the pas sage of fast (>100km s : ) shocks 



(e.g. ISanchez Contreras et al.ll2000l ; IVazquez et al.ll2000l ). From these results we infer the 
existence of fast, pAGB winds that interact hydrodynamically with the AGB CSE leading 
to the formation of shocks and, ultimately, to the acceleration and shaping of the nebular 
material. 

Direct detection of pAGB winds is limited to a very few objects. In some cases, there is 
Balmer and forbidden line e mission arising in a set of compact, shock- excited regions located 
along the nebular axis (e.g. iBujarrabal et al.lll998l ; iRiera et al.ll2006l ). These 'knots', which 
usually move away from the star at high velocity, are thought to result from the propagation 
of shocks in the pAGB wind itself, suggesting that pAGB winds are collimated and directed 
along the nebular axis. In some pPNs/PNs, fast, pAGB winds are also revealed by P-Cygni 
profiles close to the central star. One of the most exciting results recently provided by long- 
slit spectra with HST /STIS is the first direct observation of the spatio- kinematic structure 
of a fast collimated pAGB wind (jet) in a 'pristine' stage, i.e., no t strongly altered by the 
interaction with the AGB shell (ISanchez Contreras fc Sahaill200ll ). The characterization of 



the pristine jet in Hen 3-1475, which is collimated and expands at more than 2000km s" 1 at 
only £^10 16 cm from the central star, has had a direct impact on the different theories of pAGB 
wind collimation, ruling out purely hydrodynamical processes as the collimation agent in this 
object. 

Recognizing the importance of optical spectroscopy of pAGB objects we have recently 
obtained optical long-slit spectra a large sample of objects evolving off the AGB, most of 
them pPNs and young PNs, in order to characterize the pAGB mass-loss process and the 
jet-sculpting of AGB winds. This forms part of our extensive, mult i- wavelength survey 
program of imaging and spectroscopic observations of OH/IR stars ( evolved mass-losing 
stars with OH maser emission), which may be the youngest pPNs (e.g. ISahai et al.ll2007bl ; 
Sanchez Contreras fc Sahail 120041 ) . 



In this paper, we present the spectroscopic database resulting from our observations. 
Our work complements recent optical spectroscopic surveys of pPNs/PNs (e.g. ISuarez et al. 



20061 ; iPereira &: Miranda! 120071 ) by providing deep, high spectral resolution spectra covering 
a wide wavelength range. Here, we focus on the characterization and analysis of the different 
types of Ha profiles and their correlation with several stellar and envelope parameters. The 
shape of the line profile is interpreted in terms of present-day mass-loss in these objects, 
which enable obtaining information on their on-going pAGB winds. We have also derived 
spectral types based on the absorption line spectrum of our targets. Whenever available, the 
O 1 7771-5A infrared triplet have been used as a luminosity indicator. 
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2. Observations and data reduction 

2.1. Spectroscopy with ESI at the 10 m Keck II telescope 

Observations of most targets (28 sources) were carried out in 2002 and 2003, using the 
Echellete Spectrograph and Imager (ESI; Sheinis et al., 2002) in echelle mode mounted on 
the 10m W.M. Keck II telescope at Mauna Kea (Hawaii, USA). Tabled] presents a journal 
of the observations. The detector was a MIT-LL CCD with 2048x4096 squared pixels of 
15/zm. Total wavelength coverage is ~3900-10900A. The reciprocal dispersion and the pixel 
angular scale range from 0.16 to 0.30 A/pixel and from 0'.'120 to 0'. / 168, respectively, for the 
ten echellette orders (15 to 6) of ESI. The velocity dispersion has a nearly constant value 
of 11.5km s -1 pixel -1 in all orders. We used a / /5x20" slit. For extended objects with 
known morphology the slit was centered on the nucleus of the nebula and oriented along its 
main symmetry axis. For the rest of the objects the slit was oriented along the paralactic 
angle. For the object IRAS 05506+2414, we observed two different posi tions centered on th e 



compact sources labeled as Sa and Sb in the HST images presented in ISahai et al.l (l2007bl ). 
Seeing ranged between ~0'/6 and ~0'/9. Weather conditions were photometric only for the 
nights 2003, June 4 and 2004, August 11. 

All spectra were reduced using the IRAF0 package following the standard procedure. 
The long-slit spectra were bias-subtracted (using the IRAF task esibias) and flat-fielded. 
Individual exposures for each object were combined and cosmic-rays removed. One dimen- 
sional (ID) spectra were extracted using the IRAF task apall. We used extraction apertures 
large enough to include all the nebular emission (except for Hen 3-1475; see §E])- Background 
regions were defined on each side of the aperture separated by a buffer zone for the different 
orders. The background level in these regions was fitted and subtracted. Wavelength cali- 
bration was done using exposures of copper-argon (CuAr) calibration lamps and the night 
sky lines. The comparison spectra were extracted using the same aperture as the objects. 
Lines were identified in the calibration spectra and a dispersion solution was found for, and 
applied to, each of our science targets. The velocity resolution achieved (FWHM of the 
lamp lines) is ~37km s -1 for all orders. Flux calibration was obtained by observing a num- 
ber of spectrophotometric standard stars (Hiltner 600, Feige 56, LTT 1788, BD+28°4211, 
BD+33°2642, and Feige 110) using a 6"-wide slit. 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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2.2. Spectroscopy with MIKE at the 6.5 m Magellan-I telescope 

Spectra for IRAS 17150-3224 and IRAS 17440-3310 were obtained at the 6.5-m Magellan- 
I (Baade) telescope in Las Campan as Observatory (Chile ) using the Magellan Inamori Ky- 



ocera Echelle (MIKE) spectrograph ([Bernstein et aU2003l ). (We also observed IRAS 17441-2411 



with MIKE, however, the ESI spectrum for this source has a larger S/N and, therefore, MIKE 
data are not presented.) Observations were carried out on 2003, April 27-28. For the config- 
uration used in these observations, the wavelength ranges covered were ~3200-5000A and 
4900-7300 A, respectively for the blue and red arms. Each arm is equipped with a 2048x4096 
15 //m-pixels CCD. We used a l"-wide slit, leading to a resolving power of R~28,000 (22,000) 
for the blue-arm (red-arm) spectra. The images were binned 2x2 for a final reciprocal dis- 
persion of ~0.04 and 0.10 A pixel -1 , respectively for the blue and red arms, and a spatial 
scale of ~0.28 arcsec pixel -1 . Since our targets were more extended that the length of the slit 
projected onto the sky (5"), blank sky spectra with the same exposure time as for our source 
were taken right after observing each object for background subtraction. The average seeing 
was ~1" (~1"5) for the night when IRAS 17150-3224 (IRAS 17440-3310) was observed. 

For IRAS 17440-3310 and IRAS 17150-3224 the slit was centered on the nebular nu- 
cleus and the North- West lobe, respectively. The PAs along which the slits were oriented 
rotated during the exposure ±5-7° around the position quoted in Tabled] This is because 
MIKE is located in one of the Nasmyth platforms of the Magellan telescope but is not at- 
tached to the instrument rotator neither it incorporates any mechanism to compensate for 
the field rotation. This setup provides an excellent stability but implies that the position of 
the slit in the sky changes continuously as the observation progresses. 

These spectra were reduced using the MIKE pipeline developed by Daniel Kelson (see 
Kelson et al. 2000, Kelson 2003, 2007 for an extensive description of some of the numerical 
methods used). Although observations were performed under non photometric conditions, 
we corrected the data for the wavelength-dependent response of the system by acquiring 
spectra of the spectrophotometric standard stars HR 3454, HR 4468, HR 7596, HR 8634, 
and HR 7950. All the spectra were wavelength calibrated using observations of a lamp of 
Thorium- Argon taken right after observing each target. The velocity resolution achieved 
(FWHM of the lamp lines) is nearly constant for all the orders: ~12-13km s -1 for the red 
arm and ~10km s -1 for the blue arm. One dimensional spectra were extracted using the 
total length of the slit as the aperture. Our spectra probably do not include emission from 
the ends of the lobes of IRAS 17150-3224 and IRAS 17440-3310, which lied outside of the 
slit. 
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2.3. The sample 



The objects observed in this work are listed in Tables [T] and [2j In Table|2]we give the type 
of Ha profile assigned in this study (§ 13. ip together with the class of object, the spectral type 
of the central star, the fl2/f25 IRAS flux ratio, the 60/miIRAS flux (f60), the morphology 
of the optical and/or near-IR neb ula, and the chemis try. We have adopted the primary 
morphology classification system by lSahai et al.l (j2007bl ). which establishes four main classes 
of nebular shapes: bipolar (B), Multipolar (M), elongated (E), and irregular (I). We have also 
included information on two secondary structural features defined by Sahai et al., namely, 
the presence of a dark obscuring waist across the center of the nebula (denoted by "w" ) and 
direct visibility of the central star in the optical or near- infrared (denoted by "*"). Objects 
with star-like appearance in the HST images, i.e. with no nebulosity detected around them, 
are denoted as stellar (S). References for the spectral type, published high-angular resolution 
images (optical and/or near-IR), and chemistry are given in the last column of Tabled Many 
of the objects in our sa mple are listed i n the "The Torun catalogue of Galactic post-AGB 
and related objects" by lSzczerba et al.l (120071 ). where additional information and references 
can be found. 

Our sample has been mainly selected from our multi- w avelength surveys of OH/IR stars 



(e.g. iSahai et al.ll2007bl ; ISanchez Contreras Sahail 120041 ) based on their IRAS colors and 
identification of optical counterparts. OH/IR stars generally show double-peaked OH-maser 
emission lines profiles, are strong infrared sources, and are believed to be evolved mass- 
losing stars with dense CSEs. The IRAS spectral energy distributions (SEDs) of most of our 
targets show a lack of hot dust (F 12 < -F 25 ) - indicating a recent cessation of the large-scale 
AGB mass-loss process, which is believed to signal the beginning of pAGB evolution. The 
complex, aspherical morphologies of most objects as seen in high-angular resolution images 
indicate that pAGB wind interactions are taking, or have taken, place in these sources. 

Although our sample is composed in its majority of pPNs, it also includes five long- 
period variable (AGB) stars: IRAS 01037+1219 (IRC+10011), IRAS 02316+6455 (V 656 
Cas), IRAS 03507+1115 (IK Tau), IRAS 09452+1330 (IRC+10216) and IRAS 10131+3049 
(CIT 6). Their IRAS colors {Fi 2 > F 2 ^) indicate that the heavy AGB mass-loss process is 
still going on. The large-scale envelopes around these objects are roughly circularly sym- 
metric, however, the core nebular structures are often asymmetric and display bipolar or 
irregular /clumpy morphologies. This indicates that these objects are likely precursors of 
bipolar pPNs in which the jet-sculpting of the AGB CSE by interaction with collimated 
(probably fast) outflows has already started. Our sample also contains several young PNs 
(yPNs). We have adopted the classification as yPNs for objects with O- and B-type central 
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stars surrounded by newly formed, compact H II regions like the well known pAGB ob- 
jects Ml-92 and Hen 3— 1475 (see Table[2]). We have also included in our sample the object 
IRAS 19114+0002 (also known as AFGL 2343), which has a controversial evolutionary status. 
IRAS 19114+000 2 has been classified as e ither a pPN descended from low- or intermediate- 
mass progenitor (IReddy fc Hrivnakl|l999l ) or a yellow hypergiant (YHG) descended from a 



mass ive (~30M Q ) Population I supergiant (e.g. lJura fc Wernerlll999l ; ICastro-Carrizo et al. 



20071 ). YHGs are believed to follow an evolutionary path similar to that of pPNs: they 



undergo a strong mass-loss process leading to massive CSEs while their central stars become 
progressively hotter, therefore, we consider appropriate to include IRAS 19114+0002 in our 
observations. The source IRAS 05506+2414, although satisfies our selection criteria, is most 
likely not an evolved star. This enigmatic outflow source was serendipitously discovered 
as part of our multi-wavele ngth survey of pP Ns and is probably associated with an early 
stage of a massive star's life (jSahai et al.ll2008l ). The ESI spectra of IRAS 05506+2414 taken 
at position Sa (see § 12. ip is presented here for completeness but this object is not further 
discussed except for spectral type assignation (§ 14. 2p . Finally, although most of the objects 
in this study are O-rich, we have also incorporated some C-rich envelopes as well as objects 
with a mixed or uncertain chemistry. 

Our sample is not complete on any of the properties summarized in Table [2], but rather 
should be considered as a pilot sample. As happens in other studies based on optical data, 
our sample is necessarily biased toward objects with optical counterparts. This implies 
that we may be missing pAGB objects heavily obscured by thick, dusty envelopes, i.e. 
probably the most massive and youngest pAGB stars. What is the fraction of such massive, 
heavily obscured objects and what is the mass/age threshold in our sample is very difficult to 
determine because the star and envelope evolution (as well as their mutual influence) along 
the AGB and pAGB phases, which ultimately determine when the envelope becames visible 
in the optical, remain very poorly characterized. 



3. Results 

For most of the targets, our two-dimensional (2D), long-slit spectra show spatially unre- 
solved or marginally resolved emission. For the sources IRAS 05506+2414, IRAS 17423—1755, 
IRAS 19343+2926, IRAS 21282+5050, IRAS 22036+5306, IRAS 17317-2743, IRAS 17440-3310, 
IRAS 17441-2411, and IRAS 17150-3224, the nebular emission is resolved. The 2D spectra 
of these spatially resolved nebulae will be studied in detail in the future. In this paper, we 
present and discuss extracted ID spectra obtained using an aperture that includes emission 
from the whole nebula. In most cases, the spectrum is dominated by that of the bright nu- 
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cleus, the emission from the extended nebulosity being undetected or very faint. In the case 
of Hen 3-1475, the aperture has been chosen to include emission only from the bright nebular 
core (inner 2") and leaving outsi de the emission from the shocked, axial knots that are well 
separated from the nucleus (e.g. ISanchez Contreras fc Sahail l200ll ) . For the objects in our 



sample with optically thick equatorial regions (a total of 8; see references given in Table[2]), 
the light from the nucleus, including the stellar continuum, is indirectly seen scattered off 
by the dust in the lobes. 

The spectra of our targets are composed by a number of absorption and/or emission 
lines superimposed on a fairly red continuum (Fig.[T]). Most absorption lines are expected 
to be of photospheric origin whereas emission (recombination and forbidden) lines are most 
likely nebular. Other absorption features that may have an inters tellar or circumst ellar 
origin, such as the diffuse bands at 4430, 5780, 5797, and 6284A (e.g. iLuna et al.ll2008l . and 
references therein), are also observed in some of our targets. The continuum emission is 
dominated by the stellar photospheric continuum (we have checked that the contribution 
by nebular continuum is negligible, as expected, even for the objects displaying the most 
intense Ha emission). 



3.1. Hydrogen lines 

Half of the objects in our sample (15 out of 30) show Ho; mainly in emission (Fig. [2]) 
and 11 targets show Ha mainly in absorption (Fig. [3]). We have also found 3 sources (all 
of them are AGB stars) with neither Ha emission or absorption and 1 object (the pPN 
IRAS 19292+1806) in which the low S/N of the spectrum prevents us from determining 
whether Ha emission or absorption is present. 

Based on the Ha profile we can identify objects of different types (Table[2]). Among 
the sources where Ha is observed mainly in emission (Fig.[2|), we find two basic types of 
profiles, namely, symmetric and asymmetric. Only 2 sources belong to the first category 
(referred to as pure emission sources - pE): the yPN IRAS 19374+2359 and the AGB star 
IRAS 09452+1330 (IRC+10216). In CIT 6 there is strong line blending around the Ha line|J 
since Ha cannot be isolated, determining the shape of its profile is not possible. The rest 
of the Ha-emitting sources, which is also the majority (12 out of 15), show signs (in a dif- 
ferent manner and amount) of absorption bluewards of the line core, which is responsible 
for their asymmetric profiles. Most sources show clear P Cygni like profiles, i.e. with a 



2 We note that the features around Ha are not due to the noise but represent real emission lines, many of 
which are also observable in the spectrum of IRC+10216. 
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narrow/sharp blue-shifted absorption feature against the line wings. IRAS 08005— 2356 and 
IRAS 19024+0044 display two of the most remarkable P Cygni profiles amongst the objects in 
this category (referred to as pcyg). In IRAS 19520+2759 (and tentatively IRAS 17516-2525 
and IRAS 19306+1407), not one but two different blue-shifted absorption features are ob- 
served. In 3 objects (He3-1475, IRAS 21282+5050, and IRAS 20462+3416), the Ha blue 
wing is clearly weaker than the red one. Although a sharp blue-shifted absorption is not 
observed in our spectra, these objects belong to the pcyg class. In the case of the well known 
yPN Hen 3— 1475, HST /STIS optical spectra show a remarkable Ha P Cygni profile with 
two distinct blue-shifted absorption features well delineated against the line wings (Sanchez 
Contreras & Sahai 2001). These absorption features appear smoothed out in our ESI spectra 
and, more generally, in ground-based observations probably due to their blending with emis- 
sion components arising in different nebular regions superimposed within the PSF. In the 
case of IRAS 21282+5050 and IRAS 20462+3416, P Cygni like profiles are observed in other 
H I and He I lines (see e.g. H/3 in Fig.HJ, which supports Ha having a true, but partially 
masked out, P Cygni profile. Moreover, further confirmation of P Cygni line profiles exists 
from previous spectroscopic observations of both IRAS 21 282+5050 and IRAS 20462+3416 



flSmith & Lambertlll994 



Arrieta fc Torres- Peimbertll2003l . and references therein). 



In most of the objects showing Ha absorption (Fig. [3]), the line profile is found to be 
partially filled with emission (see a detailed fitting of the different absorption and emis- 
sion components of the profile in § 14.11 ) We refer to this objects as emission filled absorption 
sources (efA). The emission component is quite prominent, reaching or exceeding the contin- 
uum level, for example, in IRAS 17440-3310, IRAS 19114+0002, and IRAS 19475+3119. In 
other cases, the Ha absorption shows a peculiar, structured profile with an emerging emission 
feature or "hump" bluewards of the narrow, deep absorption core: e.g. IRAS 04296+3429, 
IRAS 17441— 2411, and IRAS 18167— 1209. Such complex emission-absorption Ha profiles 
are most likely the result of the superposition of two lines of different origin: a photospheric 
absorption core and an emission component that originates in the stellar close surroundings 
(see further discussion in § 15.31) . There are only two objects in our sample, IRAS 17150—3224 
and IRAS 19477+2401, for which the Ha line presents a pure absorption (pA) profile, at least, 
within the limited S/N of our data. In a number of efA and pA objects, the Ha absorption 
profile has two components: a deep, narrow central absorption feature and broad absorption 
wings. Similar Ha profiles are typically found in F-type supergiants, in which the broad 
wings are most likely due to the Stark broadening mechanism operating at deep layers of 
the stellar atmosphere, whereas the narrow component is consistent with being formed in 
the low surface gravity photosphere of these stars. The blue emission-hump observed in the 
absorption profile of some objects is also commonly found in F-type supergiants. 

We briefly discuss now the profile of other recombination H I lines, whenever they 
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are detected, and compare it with that of Ha. As expected, none of the 3 AGB stars 
with Ha non-detections show other H I lines either in absorption or emission. In the pPN 
IRAS 19292+1806, which is not detected in Ha because of low S/N in that spectral region, 
weak Paschen lines (from Pai6^3 at 8502A to Pa 9 ^ 3 at 9229A) are seen in absorption. For 
pA and efA objects whenever other lines of the Balmer or Paschen series are detected they 
show pure absorption profiles — i.e. the incipient emission seen in Ha for efA objects is 
not visible in other H I lines. Among the two objects for which Ha is observed with a 
pure emission profile there are no other H I lines detected except for CIT 6, which shows 
H/9 emission with a quite broad, most likely blended profile. Except for IRAS 17516—2525 
and IRAS 22574+6609 (see below), all sources with Ha P Cygni profiles display the same 
type of profile in one or more of the Balmer lines and also, in some cases, in the Paschen 
series (observable in our ESI spectra from the Paschen discontinuity at 8200A to Pa7_^3 at 
10049. 4A; note, however, that the Pa lines from upper levels n=lQ, 15, and 13 at 8502, 8545, 
and 8665A, respectively, are usually blended with the near-IR Ca II triplet and determining 
their profiles unambiguously is difficult). The P Cygni profile is more easily recognized 
and remarkable in H/3 than in Ha, in which the intense emission wings partially infill the 
absorption component (see e.g. IRAS 21282+5050; Figs f5] and H|). The absorption component 
of the P Cygni profile appears more prominent (sharper and deeper) relative to the emission 
for bluer Balmer lines. This is partially an effect of the emission component becoming 
progressively weaker for the recombination lines from higher levels (which are expected to 
be less and less populated as the quantum number n increases) . In some cases, the emission 
component almost disappears completely for the bluest Balmer lines observed by us (H7 
and H<5), resulting in a blue-shifted pure absorption profile. In contrast to the emission 
component, the absorption strongly depends on the population of level n=2, which is the 
same for all the Balmer lines. IRAS 17516-2525 and IRAS 22574+6609 are the two only 
pcyg sources for which Pa lines are observed with a pE and pA profile, respectively, i.e. non 
P Cygni. (In both cases, the S/N in the spectral region of the Balmer series is too low to 
determine whether or not these lines are present.) 



3.1.1. Variability of the Ha profile 

We have noticed differences in the Ha profiles of some of our targets with respect 
to earlier observations. Here, we describe and briefly discuss such differences for compact 
(unresolved or marginally resolved) objects in which the observed variations most likely 
represent real changes, e.g. are unlikely to be due to the different slit width and orientation 
in our spectra and earlier datasets. 
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IRAS 19114+0002. Our Ha profile shows a narrow absorption feature (A b s =6564.5A) 
and two adjacent (blue- and red-shifted) emission components (Fig.[3]). There is another 
absorption feature b luewards o f Ha ( A n hs=6561.lA) that is most likely the Ti I A6559.57A 
line as suggested by IZacs et al.l (119961 ). In fact, the center of the Ha absorption component 
and that attributed to the Ti I line yield a very similar value of the systemic velocity, 
Vlsr=1H and 104 km s _1 , respectively, in good agreement with earlier results from other 



absorption lines and C O measurements ( e.g. iReddy fc Hrivnak!ll999l ; iBujarrabal et al.ll2001 



and references therein). IZacs et al.l (119961 ) also observed a double-peaked emission component 
and a narrow absorption core in the Ha profile of this object, i.e. a shell profile, however, 
in our spectrum the red emission component is stronger than the blue one, in contrast 
to what these authors found. Moreover, the velocity difference betwee n the two emissio n 
peaks measured by us, ~145km s -1 , is larger than that observed by IZacs et al.l (119961 ). 
~128km s -1 . 

IRAS 19475+ 3119. Variations o f the Ha profile obs erved with ESI with respect to pre- 
vious measurements are reported by lSahai et al.l (j2007al ). Our Ha profile shows a broad ab- 
sor ption feature and a narr ow inverse P Cygni shaped core (Fig.[3]). The Ha profile observed 
by iKlochkova et al.l (120021 ) shows similar broad absorption wings, howev er, the line core 



shows two emission peaks neighboring a deep, narrow absorption feature. IKlochkova et al. 



(120021 ) also present and compare their Ha profiles for three different epochs and show that 
the intensity of the red (blue) emission peak gradually increases (decreases) with time. The 
line-shape observed by us fits then quite well with the observed trend, representing the 
extreme case in which the blue peak has disappeared completely. 



IRA S 2 04 62+ 34 16. Our Ha profi le is different from that previously reported by lSmith fc Lambert 



(11994J ) and iGarcia-Lario et al.l (119971 ). in particular, the blue-shifted absorption component 
of the P Cygni profile is less pronounced in our observations. These authors already no- 
ticed strong changes in the profile of Ha and other l ines in a time-scale o f days, which 
are interpreted in terms of mass- loss episodes (see also lArkhipova et al.ll200ll ). We confirm 
such changes in some recombination lines, for example, we detect the He IA6678A line in 
absorption, however, this line has had a remarkable P Cygni profile in the past. 

AGB stars: IRC+10216, CIT6, and IK Tau. Line emission in long-period variable stars 
is known to be a transient phenomenon. Ha emi ssion was not present in the spectrum of 
the C-rich variable st a r IRC +10216 observed by iTrammell et al.l (119941 ) but was detected 
by ICohen fc Schmidtl (119821 ) as well as in this work. The other C-rich AGB star in our 
sample, CIT6, has gone through emission-line phases showing intense Balmer line emission 



and other non- hydrogen lines like [O I], [N II], and [S II 



excitation shocked zone (ICohenlll980l ; ICohen &: Schmidt 



suggestive of emission from a low- 



19821 ; ITrammell et al.lll994h . In our 
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ESI spectrum, Ha appears much weaker than in earlier works and is blended with other 
(presumably metallic) lines with similar intensity. Finally, we do not detect Ha or any other 
emission lines in the O-rich Mira variable IKTau, however, this object must have shown 
line emission in the past according to its M6e-M10e assignation in the General Catalogue of 
Variable Stars (GCVS). 



3.2. Non-Hydrogen line emission spectrum 



Most objects in our sample with intense Ha emission (i.e. pcyg and pE) also exhibit 
other emission lines by heavier elements. In two cases, namely, the yPNs M 1-92 and Hen 3- 
1475, the underlying stellar absorption spectrum is totally masked out by a wealth (hundreds 
to thousands) of emission lines by different neutral and ionized atoms, like He I, N I, N II, 
C I, O I, O II, S II, Fe I, Fe II,... (see full spectra of our targets in Fig. [I]) - - a com- 
prehensive list of nebular emission lines identified in the optical sp e ctrum of M 1-92 (many 
of them also detected in Hen 3-1475) is reported by [Arrieta et al.l (120051 ). The rest of the 
Ha-emitting sources in our sample, which is also the majority, show a small or moderate 
number of emission lines, lacking many of the intense forbidden emission lines ty pical of 



evolved pNe with hot (T efT 3>30,000K) central stars (e.g. Ivan de Steene et al.lll996a Jbl). Al- 
though the spectrum of the two AGB stars in our sample with Ha emission, IRC+ 10216 
and CIT6, is clearly shaped by molecular (e.g. VO, TiO, and ZrO) absorption bands and 
metallic absorption lines typical of C-rich late type stars, several weak emission lines are also 
detected. Hydrogen and metallic emission lines are known to be prese nt in the spectrum of 



long-period variables at certain phases of their pulsational cycle (e.g. ICastelaz et al.l 12000 



and references therein). Such emission lines are interpreted as due to heating of the stellar 
atmosphere by shock waves. In these two objects we also detect very prominent C2 Swan 
bands A6059 (2,4), A6122 (1,3), and A6191 (0,2). To our knowledge, this is the first time 
that these particular bands are reported in these objects. In the spectrum of CIT6, we also 
detect the band A5097 (2,2) and tentatively A5070 (3,3), A5635 (0,1), and A5585 (1,2). The 
Swan bands of the Co molecule are comm only observed in C-rich stars and comets (see e.g. 
Klochkova et al.l Il999i ; iBieging et al.l 12006k and references therein) . 



Most sources in which Ha is observed in absorption, i.e. with a pA or efA profiles, 
have a spectrum dominated by absorption lines with no hint of nebular emission lines. The 
central stars of pA and efA objects all have F and G spectral types, therefore, the lack of 
emission lines is consistent with these relatively cool stars not being able to ionize a significant 
fraction of their circumstellar material. The are only three objects in which Ha is seen in 
absorption that also exhibit weak emission features: IRAS 04296+3429, IRAS 19114+0002, 
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and IRAS 20136+1309. (We also find tentatively one weak emission line at ~6604A in the 
spectrum of IRAS 19475+3119, which could be identified with the very low excitation line 
Ni IA6604.29A.) 

Two relatively intense emission features around 5632 and 5583A are observed in the 
compact (unresolved) nuclear region of IRAS 04296+3429 (Fig. [I]). These lines, which are 
very broad (FWZI~5A) and have a peculiar, triangular profile, are identified as the emission 
bands (0,1) and (1,2) of the Swan system of the C 2 molecule at 5635 and 5585A, respec- 
tively The (0,0) A5165, (1,1) A5129, (0,2) A6188 C 2 bands are also detected with a smaller 
S/N in our spec t rum. The (0,1) and (0,0) Swan bands have been previously reported by 
Klochkova et al.l (119991 ). whereas the other three vibrational transitions are reported by us 
for the first time. The emission band (1 , 0) A4735 is abs ent in our spectrum as well as in 
that presented by lKlochkova et al.l (119991 ). iHrivnaki (119951 ). however, report detection of this 
band in absorption from low-resolution spect ra. We c o nfirm the presence of a weak emission 
feature around 4071 A previously noticed by IHrivnaki (119951 ). This line roughly coincides in 
wavelength with an emission feature (attributed to a blend of [S II] lines) observed in the 
C-rich pPN CRL 2688. There are two other unidentified emission lines in the spectrum of 
IRAS 04296+3429 at 8065.6A and 8102.9A. The narrow profile of these emission features is 
consistent with being atomic or ionic lines rather than molecular bands. 

There are only 2 emission lines (other than Ha) observed in the spectrum of IRAS 19114+0002 
that we identify as the [Ca II] AA7291. 47,7323.891 doublet. These lines, which arise in the 
compact (unresolved) circumnuclear region, are also observed in other objects in our sam- 
ple: IRAS 17516-2525, Ml-92, Hen 3-1475, IRAS 22036+5306, and IRAS 08005-2356, and 
tentatively in IRAS 19520+2759. (We can't rule out or assert the presence of these lines in 
the spectra of IRC+10216 and CIT 6, which display many emission lines strongly blended in 
that region). Since calcium has among the highest gas phase depletions, detection of these 
lines indicate that calcium is not depleted onto dust grains either because the [Ca II] doublet 
emission arises in a dust-free environment or because Ca is liberated from the grains. We 
consider grain destruction by moderate-velocity shocks a likely explanation in the case of 
IRAS 19114+0002 and, more generally , for the pPNs and yPNs in our sample with [Ca II]- 
emission (see e.g. iHartigan et al.l 119871 ) . Moderate speed motions are in fact deduced from 
the FWHM of the [Ca II] doublet in our targets, which range between 40 and 100km s -1 . 
These values are larger than the typical expansion velocities of AGB CSEs and probably 
result from shock-acceleration produced by the pAGB-to-AGB wind interaction (§[!}). It 
is worth mentioning that the [Ca II] doublet is also obser yed in the young stellar object 
IRAS 05506, which shows a low-excitation shock spectrum (jSahai et al.l 120081 ) . The critical 
density for the [Ca II] doublet is ~10 6 cm -3 , therefore, emission of these lines arises in a 
relatively tenuous gas. 
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We detect a few weak emission lines in IRAS 20136+1309 (Fig.[T]). Some of these fea- 
tures are also present in the spectra of the yPNs Ml-92 and Hen 3-1475 and the pPNs 
IRAS 08005-2356 and IRAS 22036+5306, for example, the emission around 5955 and 7911, 
which we tentatively identify with Fe IA5956.7A and Fe IA7912.8A, respectively, or the lines 
7083 (Ti I A7084.2A or C I A7085.5A?) and 8045A ([CI IV] A8045.6A or Mg I A8047.7A?). 

4. Analysis 

4.1. Parameter analysis of the Ha-profile 

The presence of Ha emission from the compact nebular core displaying either a pE, pcyg, 
or a structured efA profile is interpreted as an indication of on-going (i.e. pAGB) mass-loss 
most likely in the form of a stellar wind (see also §|5J). In this section, we parameterize 
and analyze the observed line profiles, which is needed to derive information on the current 
stellar wind and other processes that may be affecting the observed line shape. 

4-1.1. Pcyg and pE sources 

The parameters used to describe the Ha profile for pcyg and pE sources are given in 
Table [3J The emission component of the Ha profile consists of an intense core plus weak 
broad wings. In pcyg sources, the blue emission wing is affected by absorption, therefore, in 
order to estimate the full width at half maximum (FWHM) of the emission and absorption 
components separately in these cases, we have fitted a Lorentz profile to the line emission 
core and wings. The observed Ha profile has been subtracted from the Lorentzian fit to 
retrieve the absorption line-shape. In doing so, we are implicitly assuming that the intrinsic 
emission profile as produced in the nebula nucleus (i.e. before being altered by the absorption) 
is symmetric. Except for the FWHM of the absorption and emission lines, the rest of 
the parameters in Tabled have been measured directly on the observed profile to avoid 
uncertainties resulting from the quality of the fit and/or the validity of our assumption of an 
intrinsic, symmetric emission profile. Accordingly, the equivalent widths, W\, given in the 
tabll represent the total (emission plus absorption) equivalent widths. For pcyg sources, 
W\ is then a lower limit to the equivalent width of the intrinsic nuclear emission. Given the 
typical depth of the absorption feature, the emission W\ is expected to be underestimated 



3 Although W\ is normally defined to be positive (negative) for an absorption (emission) line, we express 
both the absorption and emission W\ as a positive value for simplicity. 
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only by a small factor (less than 2). The parameters of the absorption component V and 
Vmax are, respectively, the Doppler shift of the centroid and the blue edge of the absorption 
feature relative to the emission peak, which is located at LSR V e , in km s _1 . The values 
of the full width at zero intensity level (FWZI) depend on the noise level in our spectrum 
and, therefore, they must be regarded as lower limits. Among the three objects with pE 
Ho; profiles, the presence of broad wings is only confirmed for IRAS 19374+2359. In this 
case, however, the total width of the wings is uncertain partially due to the limited S/N of 
the spectrum and the poorly characterized shape of the underlying continuum. For the pE 
source CIT 6, the Ha line is strongly blended with other emission features and, therefore, 
the line parameters, which are expected to be very unaccurate, have not been calculated. 

In order to study whether Raman scattering could explain the broad wings, we also 
fitted the profile of the Ha wings with a function of the type I a oc A~ 2 , which is expected if 
the wings are due to Raman scattering (see discussion in § 15.11) . As shown in Fig.[5j the fits 
are satisfactory for all sources except for IRAS 19306+1407 and IRAS 20462+3416, which 
display wings significantly more intense than the synthetic profile. In these two cases a 
shallower power law, e.g. of the type ocA -0 ' 7 , is needed. We note that IRAS 19306+1407 and 
IRAS 20462+3416 are among the objects with the broadest wings. 

As shown in the table, the FWHM of the Ha core emission ranges between ~50 and 
200km s^ 1 , whereas the wings reach widths of up to ±2000km s -1 . The mean velocity of the 
bulk of the gas producing the absorption ranges between V^50 and 500km s -1 , however, 
larger outflow terminal velocities of up to ~ V ma _ x ~ 800km s _1 are observed (see § 15. ip . 
No correlation has been found between the FWHM of the Ha core component and the 
FWZI or W\ of the line (Table[3]). We also investigated the correlation between the Ha-line 
parameters and the stellar spectral type but no obvious trend has been found: although it 
is true that the largest equivalent widths (W\>100A) are observed for objects with hot, O- 
and B-type central stars, very small values of W x (<20A) are found for IRAS 19306+1407, 
IRAS 20462+3416, and IRAS 21282+5050, which have central stars of similar early spectral 
types. The low values of W\ in these cases are comparable to those measured in pcyg sources 
with stars with types later than B. For these cool (A through G) stars, no detectable Ha 
emission is expected at all if photoinization were the only excitation agent in the stellar 
wind. The relatively large values of W\ found in these cases are consistent with shocks being 
an important ionization mechanism. Alternatively, the large values of W\ may indicate an 
additional source of energetic UV photons (for example a hotter companion?). In a future 
work, we will attempt determining quantitatively the importance of shock excitation in the 
stellar wind using diagnostic diagrams for other emission lines arising in the close stellar 
environment. 
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The observed values of W\ for objects with the hottest central stars in our sample are 
consistent with UV stellar radiation being the main ionization mechanism. In this case, 
the Ha-line flux is directly proportional to the number of Lyman continuum photons (Nl y ) 
emitted by the star. This is true if the Ha-emitting region is radiation bounded. Under 
this assumption, we have computed W\ for H a emission arising in a H II region a round 
the star using the ionizing fluxes provided by IVacca et al.1 (119961 ) and ISmith et al.l (120021 ) 
for different spectral types an d adopting a conve rsion factor from N Ly per second to Ha 
luminosity of 3.167xl0~ 12 erg (IBrocklehurstl Il97ll ) . The intensity of the stellar continuum 
near Ha, which is needed to calcul ate W\, has been estimated fr om the M y and R— V values 
given by IVacca et al.1 Jl996h and iMafz-Apellaniz et al.1 (120041 ) (see also |(5ox| l200oh . The 
flux t hat correspond s to th e zero i?-band magnitude, 2.19xl0 -9 ergs -1 cm _2 A _1 , is taken 
from iFukugita et al.l (119951 ). The nebular continuum has been also estimated and taken 
into account to derive W\. We find that the apparent W\ dichotomy among pcyg sources 
with hot central stars of similar spectral types (with values of Wa>100A in some cases and 
W\<2QA in others) could easily result from uncertainties in the spectral classification. This 
is because N Ly and, thus, W\ are very sentive to the stellar T eff for late-0 and early-B stars. 
In particular, for class I sources with T e g of 28.1 (BO), 26.3, and 25 kK (B0.5) we find values 
of W\ of 300, 80, and 50A, respectively, i.e. there is a factor of 6 in W\ for a difference of only 
0.5 spectral subtypes. (The steep decline of W\ with T eS is because the stellar SEDs peak 
near the Lyman discontinuity for [0,B]-type stars.) For 09 stars, we obtain predicted values 
of W A ~1000A, i.e. larger than observed (~20 and ~130A, respectively, for IRAS 21282+5050 
and IRAS 19520+2759). In the case of IRAS 21282+5050, whose central source may be a 
Wolf-Rayet [WC11] star, the discrepancy between the observed and expected value of W\ is 
sma ller since WR stars produce less Lyman continuum than non-WR stars with the same T e g 



sec 



Smith et al.ll2002l ). Values of W\ smaller than those predicted could also result in the 



case of a matter bounded H II region, for example, a non-spherically symmetric Ha-emitting 
region: if the Ha emission arises in a disk or a shell with cavities along a given axis, a large 
fraction of the ionizing UV photons will escape along the poles. 



4-1.2. EfA and pA sources 

In order to estimate the relative Doppler shift (AA), equivalent widths (W\), and FWHM 
of the absorption and emission components of the Ha line for efA+pA sources we have de- 
composed the observed profile using different Gaussian functions (Figs. [Mil and TableH]). 
Two absorption components, one narrow and one broad, are needed in most cases. For 
IRAS 19477+2401 and IRAS 19114+0002 we have only fitted one narrow absorption compo- 
nent, however, a weak broad absorption feature cannot be ruled out: for IRAS 19477+2401 
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such broad absorption may be well hidden within the low S/N of our spectrum, and for 
IRAS 19114+0002, it could be totally masked out by the relatively intense emission that fills 
in the absorption profile. For IRAS 19114+0002 and IRAS 17440—3310, the line parameters 
derived are especially uncertain, in particular, the errors in AA are probably larger than those 
quoted in the table. This is because the absorption features are almost completely filled in 
with the emission and, therefore, the fit to the different absorption components is poorly 
constrained. This leads to a poor determination of the parameters of the emission feature 
as well. Finally, the observed Ha profile for the pPNs IRAS 04296+3429, IRAS 17441-2411, 
and IRAS 18167— 1202, which shows an emission "hump" bluewards of the narrow absorp- 
tion core, can also be reasonably well fitted without the emission (line parameters with or 
without the emission component are given in TableH]). Therefore, although the line shape 
is better reproduced including a weak emission feature, the line parameters derived for this 
component are somewhat uncertain. 

Although the profile of photospheric absorption lines may be better represented by a 
Lorentz or a combined Lorentz-Gaussian (Voigt) function, we don't expect the parameters 
derived using these more complex functions to be significantly different than those in TableH] 
given the good match between the Gaussian fit and the observed line shape. Another way 
of attempting the characterization of the emission component in efA sources is to subtract 
from the observed profile a synthetic absorption line obtained from stellar atmospheric mod- 
eling. However, for an accurate characterization of the fundamental stellar parameters that 
determine the shape of the absorption profile, such as T e g, gravity, turbulence velocity, etc, 
we need to identify and fit simultaneously not only Ha but all the absorption features in our 
spectra, which is beyond the scope of this paper. 



4.2. Spectral classification: line absorption spectra 

We have used our ESI and MIKE spectra to derive the spectral type of our tar- 
gets and investigate the correlation of this fundamental stellar parameter with the type 
of Ha profile. For some of our targets, there was already an estimate of the spectral 
type available in the literature (see e.g. "The Torun catalogue of Galactic pAGB and re- 



lated objects" at http://www.ncac.torun.pl/postagb and references given therein). When- 



ever no, or uncertain, spectral classification was found, we estimated the spectral type by 
comparing the normalized spectra of our targets with those of template stars obtained 
from published stellar libraries. Note that the slope of the continuum cannot be safely 
used for spectral classification due to the uncertain flux calibration in the majority of 
the objects in our sample, which were observed under non-photometric conditions, and 
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also because these objects are expected to be very reddened by large amounts of dust 
in their envelopes. We have also double-checked the spectral types previously assigned 
to our sources. Our template spectra have been taken mainly from the Li brary of High- 



Reso lution Spectra of Stars from the UVES Paranal Observatory Project (IBagnulo et al. 



20031 . |http: / / www.sc.eso.or g/santiago/uvespop / interface.html ) . We have also used as sec- 
ondary templates some of our own targets for which a reliable spectral classification was al- 
ready assigned, e.g. from detailed studies of the stellar parameters through model-atmosphere 
methods. We have assumed luminosity class I for our objects. For several sources, this as- 
sumption is in fact confirmed from the strength of the luminosity-dependent O I infrared 
triplet (§ 14.31) . The main reference used for line identifications and laboratory wavelengths is 
the NIST Atomic Spectra Database (Version 3.1.0). In Table[2] spectral types obtained from 
this work are given using bold-face font. 

In the following we briefly discuss the spectral types assigned in this work to individual 
sources with no spectral type available from the literature and objects for which our spectral 
classification does not agree with previous assignations. 

IRAS 04296+3429. The central star o f this object wa s classified as a GO la supergiant 
based on low resolution optical spectra by iHrivnak! (119951 ). Detailed chemical analysis and 
determination of atmospheric parameters on the basis of high-resolution spectra performed 
by different authors, however, indicate a higher effective temperature of T e ff~6400-7000 K 
(IKlochkova et al.lll999l ; iDecin et al.lll998t IVan Winckel &: Reyniersll2000l ) . By comparing our 
high-resolution data with F and G UVES template standards we find that the spectrum of 
IRAS 04296+3429 is, in fact, inconsistent with a G spectral type (e.g., Paschen lines are 
extremely weak in G-type stars) but, rather, supports the classification of the central star 
as an early-F. In particular, the best match is obtained with the UVES standard HD 74180 
(F3Ia). 

IRAS 05506+ 2414- This enigmatic IRAS source, which is most likely associated with 
a YSO, was serendipitously discovered in our survey of pPNs (see § 12.31 and Sahai et al., 
2008). Its central source is listed as an M 6 in the SIMBAD database. Our ESI spectrum, 
however, is not consistent with such a late spectral type (note, e.g., the remarkable differences 
with respect to the spectra of M-type stars in our sample). The deep infrared Ca II triplet 
together with the absence of the Paschen series is indicative of a G-K star. Considering also 
the number and depth of metal lines in the 8400-8800A window, e.g. Ti I and Fe I lines, we 
assign a spectral type G9-K2 to the central source of this object. 

IRAS 17150— 3224- The central star of this object has been previo usly classified a s 
G2I based on a low-spectral resolution spectrum in the 4000-6800A range ( Hu et al.lll993 ). 
However, the density and depth of the lines across the whole wavelength range observed 
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with MIKE is most consistent with an F3-7 spectral type. (The windows 3890-4600, 5200- 
5400, 5800-6400, and 7050-7250A are particularly useful for discriminating between F and G 
spectral types.) In particular, the shape of the profile of the Ca II H&K lines (at ~3950A) 
indicates an F3-7 star. The spectrum of the template UVES star HD136537, which is a 
G2 II, has a much more dense absorption line spectrum and the lines are deeper than in 
IRAS 17150-3224. 

IRAS 17317— 2743. The central star of this object is classified as a F5 1 b ased on a low- 



resolution optical spectrum in the 4000-6900A range by ISuarez et all ( 20061 ). We observe 
very deep near- IR Ca II triplet lines (in the 8500-8700A range), which is consistent with 
an F or G spectral type. However, the absence of Paschen lines series indicates a G-type 
star. The spectrum of IRAS 17317—2743 best resembles that of the UVES template star HD 
204075, which is a G4Ib. 

IRAS 1 7440—3310. There is no previous spectral classification for this object available in 
the literature. The spectrum of IRAS 17440-3310 is quite similar to that of IRAS 17150-3224 
suggesting an F-type star. In the region of the Ca II H&K lines as well as in the 5200-5400, 
5800-6000, and 6100-6300A ranges, which are particularly useful to distinguish between 
early and late F spectral types, the spectrum of IRAS 17440—3310 is most consistent with 
an early-F spectral type. In the following we adopt an spectral type F3 1 given the good 
agreement with the spectrum of our template UVES star HD 74180 (F3Ia). 

IRAS 17516— 2525. No previous spectral type assignation is reported for this object. We 
are unable to obtain an accurate spectral type for this source since our ESI spectrum does 
not show any absorption lines but rather is dominated by nebular emission lines. The shape 
of the spectrum at long wavelengths enables ruling out an M-type classification. We note the 
remarkable similarity between the line-emission dominated spectrum of IRAS 17516—2525 
and that of IRAS 19520+2759, especially redwards of ~7400A. Since the latter is tentatively 
classified as an 09 (see below), this may suggest a hot central star also in this case. 

IRAS 18167— 1209. There is no previous spectral classification for this object available 
in the literature. Its spectrum is similar to that of IRAS 19114+0002, suggesting a mid- or 
late-F type. The shape of the profiles of the near-infrared Ca II triplet and Paschen lines in 
the range AA8400-8900 as well as the number and depth of the lines across the wavelength 
range covered by ESI (but particularly in the 5000-5600 and 6100-6300A windows) are very 
similar to those observed in the UVES template star HD 1089668, F7Iab. 

IRAS 19292+1806. No spectral classification has been previously reported. The S/N 
in our spectrum is too low in most of the wavelength range observed, including the region 
around Ha. Relatively weak Paschen lines are seen in absorption in the 8400-8900A window, 
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which may indicate a B spectral-type in this case. 

IRAS 19306+1407. Dis crepant spectral types h ave been pre viously assign e d to t he cen- 
tral star of this object by iKellv fc Hrivnakl (|2005h (B0:e) and ISuarez et all (bood ) (G5). 



Both works are based on low-resolution optical spectra. A value for the stellar tempera- 



ture o f T e gf=21,000 K, typical of a Bl type star, has also been found by lLowe fc Gledhill 



(120071 ) from detailed SED modeling leaving the stellar temperature as a free parameter. The 
absence of the near-IR Ca II triplet in our ESI spectrum argues against a G- or even an 
F-type assignation. Our ESI data support a B0-1 type central star based on the presence 
of relatively intense He I lines (e.g. at 5876 and 6678A), the density and depth of many 
metal absorption lines, e.g., in the 4500-4700 and 5600-5800A ranges, which are typical of 
early-B stars, and the relatively weakness of the Paschen lines. The presence of the C II 
lines at 6578.05 and 6582. 88A in IRAS 19306+1407 further supports its B-type classification, 
since these lines are very prominent in B-type stars but extremely weak or absent in other 
spectral types. In the following we adopt a spectral type B0-1 1 given the similarity between 
IRAS 19306+1407's spectrum and the template UVES stars HD 112272 (B0.5Ia) and HD 
148688 (Blla). Discrepant spectral types for this target (and other objects) are discussed 
in detail in AppendixlAl 

IRAS 19374+2359. The spectrum of this source is similar to that of IRAS 19306+1407, 
suggesting a B-type also in this case. The very weak Pa lines, the absence of metallic lines 
in the blue, e.g., in the AA5100-5400 region, the presence of the He I6678A line, and the very 
deep C II AA6578,6583A lines, unequivocally indicate an early B-type star. The absence of 
strong metallic lines in the 5600-5700A range and the relatively strong O I triplet suggest a 
spectral type later than B2. The spectrum of IRAS 19374+2359 best resembles that of the 
UVES template stars HD 168625 (B2/5Ia) and HD 105071 (B6Ia/Ib), therefore, we adopt 
an intermediate B3-6 1 classification. 



IRAS 194 77+ 2401. ISuarez et al.l (120061 ) classify the central star of this object as an 
F4-7I based on low-resolution optical spectra. We find, however, that the spectrum of 
IRAS 19477+2401 is rather different from those of objects with F3-F7I central stars in our 
sample (see Tabled and Fig.[T]). The spectral differences are especially noticeable in the 
8400-8900A window, in which the weak Paschen lines and strong near-IR Ca II triplet lines 
are indicative of a later spectral type. The spectrum of IRAS 19477+2401 is most consistent 
with the UVES template GO la star HD 174383. 

IRAS 19520+ 2759. There is no previous spectral classification for this object in the 
literature. Spectral typing is difficult in this case, since the spectrum of IRAS 19520+2759 
is dominated by intense nebular emission lines. The spectrum of this object shows a general 
lack of absorption lines, suggesting that it has a fairly hot star. For example, there is a clear 
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absence of lines in the wavelength window around 5700A. This region shows the presence of 
numerous lines for Bl la spectral types (see, e.g., our ESI spectra of IRAS 20462+3416), but 
not for hotter spectral types, e.g., 09 la. The line absorption spectrum of IRAS 19520+2759 
is similar to that of IRAS 21282+5050, which is classified as an 09.51. Therefore, we tenta- 
tively assign a similar spectral type, ~09, to this source. 

I RAS 20028+3910. The central star of this object has been previously classified as 



a G4 (IKelly fc Hrivnakl 12005k and references therein) . The presence of the near-IR Ca II 
triplet is consistent with both a G and F type, however, the depth of the Paschen lines 
indicates a spectral type earlier than GO. The spectrum of IRAS 20028+3910 is very similar 
to that of IRAS 04296+3429, IRAS 18167-1202, IRAS 17441-2411, which suggests a mid-F 
classification. In fact, the relative intensity of the Ca II triplet and Paschen lines is most 
consistent with an F3-7. 

I RAS 201 36+ 1 309. The central star of this object has been previously classified as a 



G0:I (IKelly fc Hrivnakl 120051 . and references therein). The spectrum of this object is similar 
to that of IRAS 20028+3910. It shows very deep near-IR Ca II triplet lines and relatively 
intense Paschen lines indicative of a mid- or late-F spectral type. The number and depth of 
many metal lines across the spectrum as well as the shape of the profile of the Ca II H&K 
lines support an F3-7 assignation. 

IRAS 22574+6609. There is no previous spectral classification for this object. Since 
IRAS 22574+6609 is very faint in the optical (V~21.3), our spectrum has a poor S/N at 
wavelengths shorter than ~6500A and no absorption features are identified in this region. 
The infrared window ~8300-8900A is, however, very useful for spectral typing: the depth of 
the Ca II triplet lines, which are very weak, and their relative intensity with respect to the 
Paschen lines constrain the spectral type to Al-6. The best-matching UVES template star 
is HD 80057, which is an Al la. 



4.3. Luminosity and distance 

The combined strength of the lines at 7771.94, 7774.17, and 7775.39A of the O IA 7773A 
infrared triplet has been established as a powerful measure of the stellar luminosity (e.g. 



Arellano Ferro et al.l 120031 . and refereces therein). These lines form as a result of upward 



transitions from the metastable 3s 5 S® level, which has a high excitation potetial of 9.15 eV, 
to the 3jjPPs 2,1 levels of neutral oxygen. The O I infrared triplet is particularly strong in 
supergiants, especially for class la objects, and is easily observable for stars with spectral 
types from middle B to early G. For later G-type stars the O I 7773A feature weakens 
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progresively (IFaraggiana et al.lll988l ). 



We have used the relationship between the absorption equivalent width of the triplet, 
W\(Q I), and the absolute visual mangitude of the star, M v , calibrated by lSlowik fc Peterson 
( 119951 ) . to derive the luminosity class of the objects for which the O IA 7773A feature appears 
in absorption. The O I triplet is partially resolved in our observations, i.e., the 7771.94 
line is observed partially separated from the 7774.17 and 7775.39, which appear blended. 
For IRAS 08005-2356, Hen3-1475, and IRAS 22036+5306, W x (0 I) cannot be accurately 
measured because the O IA 7773A feature shows a P Cygni profile and we cannot rule out 
some underlying O I emission partially filling the absorption. In addition, the O I lines may 
be blended with some other absorption features falling in that wavelength range. Therefore, 
we have not estimated My for these objects. The derived values of My (Table ED agree quite 
well (w ithin 0.2 m ) with the mean of My obtained from eqs. [1] and [2] by I Arellano Ferro et al. 



( 120031 ). (These authors obtain two separate relationships between My and the equivalent 
width of the O 17771 line, eq. [1], and the O 17774+7775 blend, eq. [2].). This agreement 
indicates that the My-W\(0 I) relationship used is valid for the domains A1-G8 and My 
from —9.5 to +0.35, i.e, for the spectral and My domain covered by most of our targets. 
One exception may be IRAS 19114+0002, with a very large value of My consistent with 
this object being a luminous, massive star. We note, however, that the relative atmospheric 
oxygen abundance of its central star has been found to be 0. 5 dex higher than that in 
normal F supergiants ( Zacs et al. 1996 : Reddy fc Hrivnak 1999h . which could result in an 
overestimate of the luminosity as derived from the My-W\(0 I) relationship. The My- 
W\(0 I) relationship may not provide reliable estimates of the luminosity either for the yPNs 
IRAS 19306+1407, IRAS 19374+2359, and IRAS 22574+6609, which have B-type central 
stars. 

The errors of My are expected to range between ~0.7 and 1.2 m , taking into account 
the errors in the coefficients of the My-W\(0 I) relationship and the error of our measure- 
ment of W\{0 I), which is expected to be <0.1A. As shown in Table[5] most objects have 
W\(0 I)>lA. These high values are restricted to Class I sources ( Thomas et al.lll979l ). im- 
plying luminosities of ~1O 4 -1O 5 L . In particular, according to these authors (see their Fig. 1) 
and taking into account the spectral type of our targets (Table^D, the large majority of our 
sources belong to luminosity class la. IRAS 20136+1309 and IRAS 19477+2401 are the two 
only objects that are most consistent with class lb. We also find that the value of W\(0 I) 
measured for IRAS 19306+1407 is comparable with, but slightly larger than, the expected 
value for a BO-lIa star. We have estimated the total luminosity (Z/boi) f° r objects in Tabled 
using the values of the bolometric correction (B.C.) that correspond to their stellar spectral 
types. 
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By comparing the position of our targets in the HR diagram (Fig. [8]) with theoretical 
evolutionary tracks for pAGB objects, we find that a significant fraction of our targets (10 out 
of 13 with W\(0 I) measurements) seem to have central stars with hig h remnant masses of 



>0.8M Q , which would point to initial masses of >3M flBloeckerlll995f ). IRAS 20462+3416, 
IRAS 20136+13090 and IRAS 19477+2401, are most consistent with lower remnant masses 
of ~O.6M , implying initial masses of ~3M . Note, however, that the large errorbars in 
the luminosity, AZog(L bo i/L Q )~0.5, together with the uncertainties arising in the evolution- 
ary models themselves and the initial-to-final mass relationship, prevent us from obtaining 
accurate values for the masses. 

We have computed the luminosity of our objects, L/D 2 , by integrating their SEDs. We 
have fitted the double-peaked SEDs by two blackbody curves, one representing the reddened 
stellar photosphere, which is mainly visible from the optical to the NIR, and the other 
representing the emission at longer wavelengths by cool circumstellar dust (Fig. [9]). The sum 
of the integrated fluxes of the two blackbodies represents L/D 2 . By comparing the values of 
L/D 2 with Lboi we have obtained upper limits to the distance, D, to our objects (TableEJ). 
The upper limits arise because the total luminosity L/D 2 derived from the observed SED is 
probably underestimated. This is mainly because our SEDs are not corrected by interstellar 
extinction, which can be quite large given the proximity of our targets to the galactic plane 
(see below). Correcting for interstellar reddening, however, is not straightforward since the 
total amount of extinction along the line of sight towards our sources depends on their 
unknown distances. Moreover, in some objects there is a clear excess of emission in the 2- 
10/im range, that is, between the star and cool dust blackbody curves, denoting the presence 
of warm dust (see e.g., IRAS 17441-2411, IRAS 19306+1407, and IRAS 22574+6609; Fig.EJ. 
Since this component is not included in our simplified two-blackbody model, the value of 
L/D 2 derived from the model SED is a lower limit to the total luminosity. 

We show now that the interstellar extinction correction factor can be significant in 
some cases, even if the total extinction is not extremely high (typical values of the inter- 
stellar extinction range between Ay=l and 4 mag). The reddening correction is particularly 
important for objects with the stellar component of the SED peaking shortwards of ~l//m. 
To illustrate this, we have corrected the SED of IRAS 20462+3416 for a value of Ay^lmag 
and adopting a A -1 power law for the extinction curve. The extinction-corrected value of 
L/D 2 is almost a factor 3 larger than that in Tabled which translates into a non-negligible 
factor of 1/1.7 for the distance. The effect of interstellar reddening is also exemplified in 
Fig. El were SEDs corrected by an average value of Ay=2.5mag are shown. 



Finally, there are additional sources of errors that may be affecting our values of D 
in Table[5] such as the typical uncertainties in the spectral classification of the central star 
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(of a few sub-types), which affects B.C. and Lbol; and the fact that the reliability of the 
M v -W\(0 I) relationship is not well tested for highly evolved objects, specially C-rich stars, 
in which the Oxygen abundance may differ significantly from solar. Also, for objects with 
an equatorially-enhanced density distribution, very common among pPNs, the luminosity 
L/D 2 computed from the integrated SED depends on the viewing angle, especially whe n a 
significant fraction of the light is scattered by the nebular dust (see e.g. ISu et al.ll200ll . for 
a quatitative estimate of viewing angle effects on the SEDs of non-spherically symmetric 
nebulae.). High polarization value s, implying mostly scattered light, are indeed observed 



for some objects in our sa mple (e.g. iTrammell et al.lll994t lGledhillll2005l ; lOppenheimer et al. 



20051 ; iBieging et al.l 120061 . and references therein) . 



4.4. Correlations 



Correlations between the type of Ha profile and other stellar and envelope parameters 
have been investigated. 

Spectral type of the central star. We have found a noticeable, and somewhat expected, 
correlation between the Ha profile and the stellar spectral type. First, there is a deficiency 
of Ha emitters among AGB stars. 3 out of 4 of the Ha non-detections are AGB stars with M 
or C spectral types. There are only two late-type stars in our sample showing Ha emission, 
namely, IRC+10216 and CIT6. In these cases, the Ha emission is quite weak [W\ <4A) 
compared with that observed in the majority of the Ha emitters (Tabled]). This weak Ha 
emission is most likely temporary, since it was absent in the spectrum of these cool objects in 
earlier observations (§ !3.1.ip . Second, as shown in Tableland Fig.[10l all objects in which Ha 
is seen in absorption, i.e. efA+pA sources, have late-type (F-G) central stars. In constrast, 
the distribution of spectral types for Ha emitters, i.e. pcyg+pE sources, is much broader, 
including objects with spectral types ranging from O to G, with a relative maximum around 
B, and also very late types such as C. Ha emission is indeed expected in objects with [0,B]- 
type central stars, whose UV radiation field ionizes the circumstellar material, producing 
an H II region observable through Balmer line emission (§ |4.1.2p . For objects with mid or 
late spectral types, lacking in photoinizing photons, a likely mechanism for producing Ha 
emission are shocks. For non-pulsating stars, such shocks probably result from the interaction 
between the inner layers of the AGB envelope, i.e. the slow AGB wind, and the fast pAGB 
winds, more recently ejected. In the two long-period variable stars with Ha emission in our 
sample, shocks could either arise in the pulsating stellar atmosphere, which is the normal 



interpretation for transient line emission in this type of objects (e.g. ICastelaz et al.l 12000 
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and references therein), or result from the AGB-to-pAGB wind interaction. We note that, in 
fact, such interaction must have started in these objects given the non-spherical morphology 
of their inner CSEs. 

NIR and IRAS colors. NIR colors are sensitive to the presence of warm (~500-1000K) 
dust. The (J-H)-(H-K) color-color diagram shows that pcyg+pE objects are redder than 
efA+pA sources (Fig. fTTj) . This is also apparent in the histograms of any of the three NIR 
colors, especially H—K and J — K also shown in the figure. We have compared the observed 
values of (J — K) with the intrinsic color (J — K) for the spectral types assigned to our 
targets and found that the larger values of (J — K) in Ha emitters do not result from 
their different distribution over spectral type with respect to efA+pA sources (Fig. [T2|) . The 
larger (J — K) color excess in pcyg+pE sources may be partially explained by a larger 
value of the total (interstellar+circumstellar) extinction. However, reddening only cannot 
explain the loci of most of these objects above the black-body line in t he (J — H)-(H — K) 



color-color diagram given the direction of the reddening vector (from ISchlegel et al.l 11998 
equivalent to a ocA -L6 extinction power law). Scattering by nebular dust could also affect the 
observed colors. However, since the scattering efficiency is larger at shorter wavelengths, the 
scattering net effect is blueing of the incident light, and that will not explain the observed 
H — K colors of pcyg+pE sources, which are redder than for a black body. While the NIR 
colors of efA+pA sources are consistent with the major contribution to the emission in the 
NIR bands being the reddened stellar photosphere, the position of pcyg+pE sources in the 
(J — H)-(H — K) color-color diagram highlights the presence of warm dust. This warm dust 
component is in fact confirmed in some cases by a clear emission excess between 2 and 10/iin, 
whenever such data exists (see SEDs of, e.g. IRAS 19306+1407 and IRAS 22574+6609 in 
Fig. [9]). The warm dust component must be near the star and, therefore, has probably been 
formed recently. Then, the presence of warm dust could be an independent confirmation of 
substantial present-day mass-loss in pcyg+pE targets, which is evidenced in the first place 
by their intense, circumnuclear Ha emission. Alternatively, the warm dust component in 
pcyg+pE sources could be located in a long-lived disk around the central source rather 
than in the pAGB wind itself. Finally, we note that the NIR colors of the efA objects 
IRAS 17150-3224, IRAS 17441-2411, and IRAS 20136+1309 also suggests dust grains with 
higher temperatures than those in their detached, much cooler, AGB CSEs (see also full 
SEDs for the two later in Fig. [9]). 

We have compared the location of our objects with the larger, un-biased sample used 



by iGarcia-Lario et al.l (119971 ). Our pcyg sources fall in regions III and IV defined by these 
authors, which suggest that the NIR colors of objects in such regions may result from recent 
mass loss leading to a hot dust component around the star. For a few of their objects for 
which optical spectroscopy was available, P-cygni Ha profiles were found. Our work has 
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demonstrated that pcyg sources systematically lie above the black-body line in the NIR 
two-color diagram, displaying colors that are consistent with the presence of hot dust in 
the stellar vicinity most likely resulted from pAGB mass-loss. There are two pcyg objects, 
IRAS 20462+3416 and IRAS 19306+1407, with black body-like NIR colors (Fig.rjTj) that can 
be explained just by reddening of the stellar atmosphere. For these two pcyg objects, there- 
fore, there are no indications of the presence of hot dust in the stellar vicinity although their 
P-cygni Ha profile indicates current mass-loss. In IRAS 20462+3416, significant spectral 
variations with time in th e Ha line have been report ed, pointing to episodic mass-loss events 



I3.1.ip . As suggested by iGarcfa-Lario et al.l (119971 ) the absence of significant NIR excess in 



this case could indicate that the duration and intensity of such mass-loss episodes are not 
enough to efficiently form and maintain hot dust grains in the circumstellar envelope. 

We have also found a correlation between the type of Ha profile and the IRAS colors, 
which are most sensitive to the cool (~100-200K) dust component that is remnant of the 
mass loss on the AGB (Fig. fTB"]) . The IRAS colors represented in this figure are defined as 
[12]-[25]=2.5Zo#(F 25 /F 12 ) and [25]-[60]=2.5Zo#(F 60 /F 25 ) with F flux density in Jy at the 
12, 25, and 60/im IRAS filters. As can be seen, efA+pA sources display slightly larger values 
of [12] — [25], i.e. the slope of the SED between 12 and 25/im is steeper, than for pcyg+pE. 
This is consistent with a well detached envelope lacking significant amounts of warm dust in 
most efA+pA sources, which is indicated as well by their NIR colors (see above). The two 
groups of Ha profile sources clearly segregate in the two-color IRAS diagram: the majority 
of the pcyg+pE sources (~64%) lie above the black-body line, whereas most efA+pA targets 
(~73%) are below. This separation suggests different properties of the circumstellar dust 
in the AGB CSE among the two groups. The distribution of the sources below the black 
body line, which means that the emission at the longer wavelengths is smaller than expected 
for a black-body, can be explained by the dependence on wavelength of the dust optical 
depth, r v oc v& (normally /3~l-2) and the fact that the dust emission is proportional to 
B u (T d )x(l — e _r "), with B u (T d ) being the Plank function for a given dust temperature T d . 
For optically thin dust emission the flux is, therefore, proportional to B v (T d ) xz/' 3 (this is the 
so called "modified" black-body function), which explains the smaller emission at the longer 
wavelengths (Fig. fTBl . The location of the targets above the black body line in the IRAS 
color-color diagram means, as for the NIR, that their colors cannot be well represented by 
a black body or modified black body with a single temperature. Such location must result 
from a larger temperature distribution of the dust over the envelope in pcyg+pE than in 
efA+pA sources. In particular, the observed effect suggests a significant contribution to 
the total emission by cold dust, especially at longer wavelengths (60/xm), whereas most of 
the emission at shorter wavelengths (12/xm) would be produced by a different component of 
warmer dust grains. 



-27- 



Galactic coordinates. We have not found any obvious correlation between the Ha profile 
and the galactic coordinates of our targets (Fig. lHj) . A correlation with the galactic latitude, 
b, would have been proof of differences in the masses of the progenitor of both types of objects. 
We note, however, that most of our targets happen to be concentrated in a relatively small 
range of latitudes around 6=0°, so the dependence of Ha with b cannot be properly studied. 
(The narrow range in b in our sample results from the fact that OH/IR surveys, from which 
most our targets were selected, normally concentrate around the galactic plane.) 

Nebular morphology. We find a variety of nebular morphologies in our sample but no 
obvious segregation with respect to the Ha profile has been found. There is nearly the same 
fraction of E, M, B, I, and stellar sources as well as w- and * objects among Ha emitters 
and efA+pA targets (see Tableland references therein). 

Chemistry. It seems that there is little connection, if any, between the chemistry and 
the characteristic of the Ha profile. There are more or less the same fraction of O- and C-rich 
objects among the Ha emitters and efA+pA sources. Among AGB stars, we find that the 
only objects with Ha emission are C-rich. Because of the small number of AGB stars in 
our sample and the expected variability of the Ha emission with the light curve, this trend 
cannot be regarded robust result. 

5. Discussion 
5.1. Fast pAGB winds probed by P Cygni like profiles 



P Cygni profiles are known to be indicative of on-going mass-loss in the form of a 
stellar wind. The formation of typical P Cygni profiles can be understood qualitatively by 
a simple model of a spherically symmetric wind in which the velocity increases outwards up 
to a terminal velocity of the outflow v^. The wind material in front of the star is moving 
towards the observer and, therefore, absorbs the stellar continuum producing a blue-shifted 
absorption feature with a Doppler shift between — and 0, relative to the source systemic 
velocity (K ys ). Except for the part that is occulted by the star, the rest of the wind, which 
forms a halo around the star, will produce broad line emission centered around V sys with 
wings extending up to iiw The sum of the emission and the absorption at each Doppler 
velocity range leads to a so called P Cygni profile. These P Cygni profiles are often seen in 
UV resonance lines of the central stars of PNs (e.g. |Perinottolll993l ; iKudritzki et al.l 119971 ). 



The Ha profiles observed in the majority of our targets differ from the typical P Cygni 
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profiles described above. The emission component is not centered at V^ ys but, rather is red- 
shifted in all cases, except for IRAS 21282+5050 (Table [3]). Moreover, the terminal velocity 
of the outflow, as measured from the blue-shifted edge of the absorption feature, is smaller 
than that obtained from the edge of the red emission wing, and in some cases, e.g. the 
yPN Ml-92 (see below), the blue absorption is well detached from the line emission core. 
These peculiar Ha line-shapes can be explained if the emission and absorption form in two 
distinct nebular components, with different velocity fields and, more generally, affected by 
different broadening mechanism. In particular, we propose the following model: the broad 
line emission, characterized by a given intrinsic (not necessarily P Cygni) profile, arises from 
a compact central source; the line emission (and stellar continuum) is scattered by dust in 
the walls of the nebular lobes, which should produce an overall red-shift of the line profile as 
a whole relative to V^ ys as observed. The blue-shifted absorption is due to neutral or partially 
ionized outflowing gas inside the lobes absorbing the scattered photons along the line of sight. 
This results in a P Cygni like profile where the blue-shifted absorption is produced against 
the intrinsic emission profile. Such a model has been successful in producing a detailed fit to 
the spatio-kinematic distribution of the blue-shifted absorption and the red-shifted scattered 
line core in the Ha line profiles observed with STIS/HST in the yPN Hen 3-1475 (Sanchez 
Contreras & Sahai 2001; see their Figure Id f or a sketch of our geometrical model). A similar 



scenario is proposed by lArrieta et al.l (120051 ) to explain the P Cygni like profile of Ha and 
other recombination lines in the yPN M 1-92, which shows a very broad, blue absorption 
feature well detached from the emission component that cannot result from the single stellar 
wind scenario used to explain typical P Cygni profiles. Our interpretation is also consistent 
with the fact that most pPNs exhibit intrinsic polarization of the con tinuum as well as som e 



emission lines, which are most likely due to scattering off dust grains (iTrammell et al.lll994l ). 
In particular, Balmer lines seem to come from an H II region close to the star, which is 
blocked from direct view by, e.g., the thick equatorial dust waists c ommonly present in these 



objec ts, and are seen in reflection off the nebular dust (see also, e.g. JSanchez Contreras et al. 



20021 ). 



The nature of the outflow producing the absorption in pcyg sources is unknown, however, 
as in the case of Hen 3-1475, we argue that it is most likely a "pristine" pAGB wind from (or 
near) the central star that has not been strongly altered by its interaction with the progenitor 
AGB wind. The projected speed of the bulk of the material in the pAGB outflow, derived 
from the centroid of the blue-shifted absorption feature relative to the systemic velocity in 
our sample, is approximately ~100km s -1 (Table[3]); however the absorption features are 
quite wide and reach larger velocities at the edge indicating the presence of smaller amounts 
of material expanding as fast as, e.g., ~800km s _1 in the case of IRAS 19520+2759. As 
for Hen 3-1475, it is plausible that a radial and or latitudinal velocity gradient exists in 



-29- 



the outflows and/or that distinct wind components are present; the latter is suggested by 
the two different absorption features observed in the Ha profiles of a few objects (§ 13. ip . 
Unfortunately, ground-based observations do not allow us to conclude on the isotropic or 
collimated nature of the pAGB winds probed by the observed P Cygni profiles, however, 
it is plausible that these are bipolar/multipolar ejections, as in Hen 3-1475 and, probably, 
Ml-92. 



5.2. Broad Ha emission wings 



The emission component observed in the Ha profile in our pcyg sources consists of 
an intense core plus weaker broad wings that reach widths (measured as the line FWZI) 
of up to ~4000km s -1 . The yet unclear origin of the v ery extended Ha wings, which ar e 
also observed in a number of PNs, has been discussed by lArrieta fc Torres- Peimbertl (120031 ). 
These authors consider several possible mechanisms for line broadening and conclude that 
Raman scattering of photons with wavelength close to Lyj3 by neutral hydrogen is the most 
probable one in 12 of the 13 objects in their sample. In this scenario, Ly/3 photons with 
a given veloc ity width are converted to optica l photons that fill the Ha wing region. The 
conclusion of lArrieta &: Torres- Peimbertl (120031 ) is mainly supported by a) the fit of the Ha- 
wing profile to a oc A -2 law; b) the presence of other Raman features produced in emission; 
and c) the presence of a significant neutral hydrogen component. Criteria a) to c) are to be 
satisfied for the line wings to be attributable to Raman scattering. 

We have checked whether or not the same conclusion can be drawn from our sample. As 
we have shown in Section I4~T| the Ha wings observed follow a oc A~ 2 law in all cases, except 
for IRAS 20462+3416 and IRAS 19306+1407. In these two objects, which happen to display 
Ha wings with similar widt hs and shapes, Raman scat t ering can be ruled out as the main 
line broadening mechanism. lArrieta &: Torres- Peimbertl (120031 ) found the same disagreement 
between the observed and predicted Raman profile for IRAS 20462+3416, which was also in 
their sample, concluding that a strong stellar wind is the most probable mechanism for 
line broadening in this case. We have searched for other emission features produced by 
Raman scattering in our spectra. In particular, we looked at the 6830 and 7088A bands 
that correspond to the Raman-scattered O VI AA1032,1038 doublet. Only four objects, 
IRAS 08005-2356, Ml-92, Hen -1475, and IRAS 22036+5306, show emission around 6830A 
with FWZI~150km s -1 , and only in two of them, Hen -1475 and IRAS 22036+5306, a weak 
emission feature is tentatively detec t ed ar ound 7088A with FWZI<100km s _1 (Fig.[T]). In 
contrast, lArrieta fc Torres-Peimbertl (120031 ) find these and other Raman features with large 



FWZIs of a fewxl0 2 -10 3 km s 1 , i.e. comparable to the width of the Ha line wings, in most 
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of their objects. 

Since both the A~ 2 profile of the wings and the presence of a large column density of 
neutral hydrogen are necessary but not sufficient conditions for attributing the observed 
Ha-wings to Raman scattering, we cannot conclude that in our sample this is the main 
mechanism for line broadening except maybe for the four objects listed above showing other 
Raman-scattered emission features. Note, however, that the FWZIs of such features are 
significantly smaller than those measured in Arrieta's sample and smaller than the Ha wings 
themselves. 



As suggested by lArrieta &: Torres- Peimbertl (120031 ). further investigation of the effect of 



Raman scattering producing the line wings could be obtained by measuring the widths of 
UV lines produced in the same region as the Ly/3 photons. This is because the width of 
the scattered Ha is proportional to the initial width of the Ly/3 line, the Raman-scattering 
broadening factor being of Xu a / AL y /3=6.4. We searched the IUE archive for UV spectra of 
our targets around the line Si III]A1892A, which presumably arises in the same region as Ly/3. 
Unfortunately, such data only exist for four objects and only in one case, Hen 3-1475, the Si 
III] line is clearly detected. The line, which was observed with the low spectral resolution 
(AA~6A) SWP, is found to be spectrally unresolved. Nevertheless, for Raman scattering 
being able to produce the observed ~4000km s _1 -wide Ha wings in Hen 3-1475, the Ly/3- 
emitting region must be characterized by a typical velocity of ~600km s -1 , i.e. the Si III] 
line should be spectrally resolved with a FWHM of ~7A as observed with SWP, in constrast 
to what is found. Therefore it seems that Raman scattering is not the main mechanism 
affecting the profile of the Ha wings in this case. 

One additional requirement for Raman scattering to be efficient is that a relatively 
strong incident Ly/3 flux has to be produced by the central source. This may be satisfied 
by Arrieta's sample, formed by objects with central stars with T eff >20,000 K, except maybe 
for the hypergiant IRC+10420 with a stellar spectral type of A5 (T e fj~8500K) at the time 
Arrieta's data were obtained]^ An estimate of the Ly/3 luminosity expected around the 
central star of the PN IC4997 using the photoinization code CLOUDY has been done by 
Lee & Hyung (2000). In this case, for a very hot (T eff =60,000 K) and quite small (R*=R Q ) 
central star, a compact H II region of high density (sil0 9 -10 10 cm -3 ) together with a column 
density of neutral hydrogen of 7VHi=fewxl0 20 cm -2 in the scattering region are needed in 
order to have a sufficient number of Ly/3 photons and efficient Raman scattering. In our 
sample, a significant number of objects showing broad Ha wings have relatively cool central 



4 The spectral type o f the central star of IRC+10420 changed from F8 to A5 in 4 years, from 1992 to 1996 
(jKlochkova et al.lll997t ). 
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stars, with spectral types F-G (T c fj~5000-7000 K) and, therefore, are unlikely to produce 
enough UV photons around Ly/3. This is in fact supported by their featureless nebular 
spectra, which are remarkably different from those observed in more evolved PNs (§ 13.21) . For 
the objects in our sample with the hottest central stars (with late-0 and early-B spectral 
types, T e fj~20, 000-30, 000 K) estimates of the Ly/5 flux will be presented in future individual 
papers. We will also investigate whether or not shocks resulting from the wind interaction 
process are able to provide a strong Ly(3 flux. 

The presence of fast stellar winds is another plausible mechanism to explain th e broad 



Ha wings observed. As already pointed out by lArrieta fc Torres- Peimbertl (120031 ). this is 
indeed the most likely mechanism for IRAS 20462+3416, for which the Ha wing profile was 
found not to follow a oc A -2 law. For similar reasons, a fast (~1400km s _1 ) stellar wind 
is also very likely being expelled by the central star of IRAS 19306+1407 leading to its Ha 
wings. In more than one third of the objects in our sample, the presence of stellar winds 
is, in fact, confirmed by P Cygni profiles in Ha and other lines. The question is whether 
or not these winds reach much higher velocities (of 1000-2000km s _1 ) in certain regions 
that can explain the broad wings observed. In the case of M 1-92 and Hen 3-1475 there 
is direct evid ence for fast winds with expansion velocities of up to ~800 and 2300km s -1 , 



respectively (lArrieta et al.l 120051 ; ISanchez Contreras Sahail l200ll ). so it is plausible that 



these winds are also responsible (at least partially) for their Ha wings. In the rest of the 
objects in our sample, the terminal speed of the fast wind measured from the blue-shifted 
absorption feature of the Ha P Cygni profile (v^ ~ V^n ax ~10 2 km s -1 , Tabled]) is never 
as large as that implied by the broad wings if these were to be attributed to fast outflows. 
However, we find that Vm ax is normally larger than the FWHM of the Ha core indicating that 
indeed there are winds moving faster than the region emitting the bulk of the Ha line (with 
an average velocity of ~120km s -1 ; Table[3]). This clearly suggests the presence of velocity 
gradients in the stellar winds of these objects. Moreover, we have found a weak correlation 
between and the FWZI of the broad Ha wings: the sources with the largest values 
of Kmx, namely, IRAS 19520+2759, IRAS 20462+3416, Ml-92, and Hen 3-1475, exhibit the 
broadest wings, reaching widths of ~2600, ~2800, ~3000, and ~4000km s _1 , respectively. 
(For IRAS 20462+3416 and Hen3-1475, the Ha absorption feature is partially masked out 
by emission, therefore, values of V mauX =— 300km s _1 are measured from their H/3 P Cygni 
profiles; Fig.Hl) This could be an indication that broad Ha wings are partially due to fast 
stellar winds. 
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5.3. Incipient mass-loss probed by efA profiles 



Structured Ha absorption profiles partially filled in with emission are co mmon in F- and 



G-type supergiants and are considered a s a possible evidence for mass-loss (e.g. iTamura fc Takeuti 



19931 ; lArellano-Ferrol Il985l ; ISowelll I1990T ) . Although most efA sources in our sample are sur- 
rounded by extended nebulosities (see Table[2] and references therein) , the observed struc- 
tured Ha profile, in particular, the weak emission that fills in the stellar absorption, arises 
in the bright nebular nucleus, i.e. in a compact region very close to the central star, and 
not in the extended nebulae. This compact nuclear region is directly seen but is spatially 
unresolved in our ESI spectra for IRAS 04296+3424, IRAS 17440-3310, IRAS 18167-1209, 
IRAS 19114+0002, IRAS 19475+3119, and IRAS 20136+1309. Among these, IRAS 17440-3310 
is the only efA source for which our ESI spectrum also shows weak Ha emission arising in the 
extended lobes (we note, however, that the spectrum of this object is yet dominated by the 
bright nucleus). In this case, therefore, part of the emission filling in the absorption profile 
in Fig.[6]may be locally produced in the lobes. Some contribution by nebular Ha emission 
to the efA profile cann ot be ruled out ei ther for IRAS 04296+3429, which is known to have a 



~l'/5 extended nebula (jUeta et al.ll2000l ) that is unresolved in our long-slit spectrum. For ob- 
jects with optically thick equatorial regions, the nuclear region is blocked from direct view. 
This is the case of IRAS 17317-2743, IRAS 17441-2411 and IRAS 20028+3910. Our 2D 
spectra of IRAS 17441—2411 and IRAS 17317—2743 show, in fact, continuum emission from 
two extended lobes separated by a dark equatorial region but there is no Ha emission locally 
produced in the lobes. The similar spatial distribution along the lobes of the stellar contin- 
uum and the structured efA Ha profile indicate that both the continuum and the Ha line 
originate at the nebular core and are scattered by dust in the lobes. This is most likely the 
case of IRAS 20028+3910, for which our ESI spectrum shows unresolved continuum emission 
arising in the brightest of its two reflection lobes. 

As shown in Section l4TTj the emission filling in the absorption in efA sources is well 
reproduced by a Gaussian function with an average width of FWHM~65km s _1 (in deriving 
this value we have ignored objects for which the width of the emission component is uncertain; 
see TableH]). After deconvolution with the spectral resolution in our data (37km s -1 , §[2]) 
and thermal line broadening (~15km s" 1 for a kinetic temperature of 10 4 K), this value 
implies moderate speed (~50km s _1 ) motions in the Ha-emitting region. This value is larger 
than the typical expansion velocity of the slow (5-15km s _1 ) AGB wind but smaller than 
the velocity of the fast pAGB outflows observed in pcyg sources (typically >100km s _1 ; 
see § 15.11) . For the majority of the objects with a reliable estimate of the Doppler shift 
(AA) for the emission feature (§ 14.11) . the emission is displaced bluewards from the narrow 
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absorption core: the blue-shift found ranges from —10 to —43km s _1 (TableHJ). Only in one 
case, IRAS 20136+1309, the emission and the narrow absorption components are centered 
around the same wavelength within the expected errors in our measurements. Among the 
three objects with the weakest Ha emission features and, therefore, most uncertain emission 
line parameters, two show blue-shifted emission and one, IRAS 04296+3429, red-shifted (by 
+7km s -1 ) emission. We note, however, that for the latter, the parameters derived from the 
Gaussian analysis of the Ha line are somewhat uncertain due to the weakness of the emission 
component and the fact that some contamination of the efA profile by nebular Ha emission 
cannot be ruled out. 

We interpret the weak emission partially filling the efA Ha profile, which arises in the 
stellar vicinity, as an indication of present-day (i.e. pAGB) mass-loss. (The material in the 
extended nebulosities around most objects resulted predominantly from the large-scale mass- 
loss process during the AGB - §[TJ) For some of our efA sources, the presence of nebular 
material near the star is confirmed by detection of a small number of weak non-Hydrogen 
emission features, e.g. IRAS 04296+3429, IRAS 19114+0002, and IRAS 20136+1309 (and, 
tentatively, IRAS 19475+3119; § 13. 2ft . Since the stellar spectral types of the objects with 
efA Ha-profiles are F and G (i.e. their central stars are relatively cool, with temperatures 
ranging between 5000-7000 K), the emission component is likely to be formed in the de- 
excitation region behind a shock wave, presumably produced by wind interaction, rather 
than in a photoionized H II region. The fact that the emission is systematically doppler 
shifted bluewards with respect to the narrow absorption core can be explained if there is 
significant occultation of the receding part of the outflowing stellar wind. Note that for stellar 
occultation effects to be significant the size of the Ha-emitting region has to be comparable 
to the radius of the star. 

Alternatively, the weak Ha emission observed in efA sources could arise, in principle, in a 
long-lived reservoir of gas near the star, such as a rotating disk, rather than in a stellar wind. 
We believe, however, that the systematic blue-shift of the emission would be more difficult 
to explain under this scenario since a long-lived circumstellar structure would necessarily 
have to be rotating (to persist throughout a substantial amount of time) and, in this case, 
there is no obvious reason why the emission from the receding edge of such structure would 
be preferentially depressed/reduced. 

Finally, although, both the narrow and broad absorption components of the efA and 
pA Ha profiles have most likely a photospheric origin, they are expected to arise in differ- 
ent layers of the stellar atmosphere, the broad feature probably arising in deeper, denser 
regions. The observed Doppler shifts between both absorption components, therefore, could 
result from the complex kinematics across the stellar photosphere. Complex, large ampli- 
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tude atmospheric motions, involving e.g. the presence of shock waves propagating across the 
atmosphere, are derived from asymmetrical absorption lines and rad ial velocity measure- 



ment s in supergiants, including those that are pPN candidates (e.g. iLebre &: Gilletl [1991 



19921 ). Such complex motions are also evidenced by variable Doppler shifts of the Stark 



photospheric absorption relative to V^ ys (see references above). 



6. Differences between pcyg and efA sources 

Both pcyg and efA Ha profiles are interpreted in terms of on-going pAGB mass loss. 
This stresses the fact that mass-loss does not stop after the AGB and that theoretical evolu- 
tionary models need to include pAGB mass loss in their models. The main difference between 
the Ha emission in pcyg and efA sources is that the strength of the line is much larger in the 
former (TableEJ&HJ). Interpreting the strength of Ha in terms of mass-loss rate is extremely 
difficult because the former depends critically on several parameters that are unknown, such 
as the relative contribution by UV photons and shocks to the ionization of the Ha-emitting 
region, whether such a region is radiation or matter bounded, and the density distribution, 
geometry, and velocity field of the gas in the stellar vicinity. Also, accurate modeling of the 
stellar atmosphere would be needed to characterize and substract the contribution by the 
underlying absorption profile to the observed Ha line shape - this is particularly significant 
for objects with A- and B-type central stars. A simple, perhaps naive, interpretation is that 
the more intense Ha emission in pcyg sources implies larger pAGB mass-loss rates. A larger 
amount of circumstellar material resulting from larger pAGB mass-loss rates in pcyg sources 
would naturally explain their NIR colors, which independently point to the presence of warm 
dust in the vicinity of the star most likely produced by substantial present-day mass loss. 

In principle, one could think that larger pAGB mass-loss rates would imply larger pro- 
genitor masses and, therefore, larger luminosities in pcyg sources. The location of our targets 
in the HR diagram (Fig. [8]) does not show particularly high values of the luminosity for pcyg 
objects. The similar distribution of pcyg and efA targets in the Galaxy does not reveal 
either a significant difference in the mass of their progenitor, although a larger sample cov- 
ering a broader range in galactic latitude would be needed to attain definitive conclusions. 
If confirmed, a lack of mass/luminosity segregation between the two types of sources would 
indicate that the unknown mechanism that governs mass-loss during the pAGB phase is not 
strongly dependent on intrinsic stellar characteristics but rather it may depend on extrinsic 
properties like, for example, the presence of a binary companion. 

Another difference between pcyg and efA sources concerns to the velocity of their current 
stellar winds (§ 15.31 and 15.11) . For efA, the weak Ha emission profile indicates moderate 
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speed motions in the nuclear Ha-emitting region, with V^ xp ~50km s _1 . For pcyg sources, 
the values of the FWHM of the Ha emission core implies larger expansion velocities in 
the stellar vicinity, typically K xp ~120km s _1 . Most importantly, the Doppler shift of the 
absorption in the P Cygni profile suggests very fast ejections reaching velocities of up to 
fewxl0 2 -10 3 km s _1 . This is consistent with the earlier spectral types of pcyg sources since 
the escape velocity (and, hence, the stellar wind velocity) is expected to increase as the star 
evolves from the AGB to the PN phase due to the shrinking stellar radius. The faster winds 
and earlier spectral types of pcyg sources could either represent a difference in age or in the 
speed traversing the pAGB path with respect to efA objects: pcyg sources seem to be older 
or to have evolved faster than efA. A more rapid evolution could again suggest more massive 
progenitor for pcyg objects, but, as mentioned above, no signs of mass variance are found in 
our sample. A faster evolution could also result from larger pAGB mass loss rates in pcyg 
sources, which would be consistent their strong Ha emission. 

With regard to the extended optical nebulosities, which are detected for most sources in 
our sample (Tabled, both classes of objects show many signs that the jet-sculpting process 
responsible for the shaping of the AGB envelope to its current aspherical morphology is 
active or has been active in the past. However, we have not found any indications of pcyg 
sources being older (nebula-wise) than efA targets since there are no obvious differences in 
their nebular morphologies (e.g., in the diversity of shapes and/or large- and small-scale 
structural components, direct visibility of the central star, or angular size). Nevertheless, for 
a proper evaluation of nebular aging, reliable estimates of the distance as well as accurate 
characterization of the nebular velocity field are needed. Finally, we would like to note that 
the on-going pAGB ejections probed by the nuclear Ha emission in pcyg and efA sources 
are probably not the same that shaped and accelerated the much more extended nebular 
lobes. We don't know whether pAGB winds (similar or not to the current ones) have been 
ejected in a continuous or episodic manner since the shaping of the AGB envelope began. 
The fact that pcyg and efA objects with central stars of similar late types (F and G) exhibit 
Ha emission profiles that differ both in intensity and shape may suggest episodic pAGB 
ejections: pcyg would represent targets that are caught at times of intensive jet activity 
leading to energetic shocks and intense Ha emission, whereas efA could be currently in a 
status of "mild" wind interaction (resulted from smaller velocities and/or densities in the 
pAGB wind?). Variable pAGB stellar winds are indeed evidenced by temporal changes in 
the intensity and profile of some nuclear emission lines, including Ha, in certain objects 
in our sample, namely, IRAS 19475+3119, IRAS 19114+0002, and IRAS 20462+3416 (see 
§ 13.1.11 and references therein). In pA targets, the pAGB winds that presumably shaped 
their aspherical AGB envelopes may have temporarily or permanently ended. 

Apart from differences in the on-going pAGB ejections, the divergence of pcyg+pE and 
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efA+pA sources in the IRAS color-color diagram suggests dissimilarities in the properties of 
the detached dust envelope produced by the large-scale mass-loss process during the AGB. 
The separation of both types of Ha profile sources in the IRAS color-color diagram seems 
to indicate a larger range of dust temperatures for pcyg objects. The different range in 
dust temperatures may result from: 1) a different distribution of the grain density in the 
AGB CSE; 2) a different spatial distribution of grain sizes; 3) a different grain composition, 
which will result in a different power-law dependence of the absorption/emissivity efficiency 
with wavelength; or any combination of these situations. For example, large temperature 
gradients are expected for optically thick envelopes since the outside grains will remain cooler 
than for optically thin CSEs due to the fact that optical-UV radiation is heavily attenuated 
and therefore unable to heat up the outside grains. Envelopes with an asymmetric dust 
distribution with much larger optical depth along one preferred direction (e.g. along the 
equatorial plane) could also result in a broader distribution of dust temperatures since grains 
in regions strongly shielded from the star radiation will remain cooler than those that are 
more exposed. Detailed modeling of the SED, including mm- and submm-wavelength data 
as well as IR spectroscopy, would definitively help in improving the characterization of the 
dust envelope (including its chemistry) and mass-loss history of these objects throughout 
the AGB and pAGB phases. 



7. Summary 

We present echelle long-slit optical spectra of a sample of evolved intermediate-mass 
stars in different evolutionary stages: 5 AGB stars, 17 pPNs, and 6 yPNs. Our sample also 
includes the object IRAS 19114+0002, which has a controversial classification as a pPN or 
a yellow hypergiant, and one YSO, IRAS 05506+2414, which was serendipitously discovered 
in our multi-wavelength survey of pPNs. (The spectrum of the latter is presented here for 
completeness but it is not discussed except for spectral typing of its central star.) We have 
analyzed extracted ID spectra of our targets with special focus on the characteristics of the 
Ha line profile arising in the vicinity of the central source, i.e. the nebular nucleus. In this 
section, we summarize the main results obtained from this work: 

- Fifteen objects in our sample show relatively intense Ha emission, whereas eleven 
targets show Ha mainly in absorption. We have also found three sources (the AGB 
stars IRC+10011, V656 Cas, and IKTau) with neither Ha emission or absorption, and 
one object (the pPN IRAS 19292+1806) in which the low S/N of the spectrum prevents 
us from determining whether Ha emission or absorption is present. 

- Based on the shape of the Ha line, we have defined four main types of sources. Among 
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the Ha emitters, those with a symmetric Ha profile are referred to as pure emission 
sources (pE), whereas objects with asymmetric P Cygni like profiles are denoted as 
pcyg. Objects with a pure absorption Ha profile are named pure absorption targets 
(pA). Objects with an absorption profile partially filled with weak emission are referred 
to as emission filled absorption sources (efA). 

The presence of Ha emission from the compact nebular core displaying either a pE, 
pcyg, or an efA profile is interpreted as an indication of on-going (i.e. pAGB) mass- 
loss most likely in the form of a stellar wind. The observed Ha profiles have been 
parameterized and analyzed to derive information on the current stellar wind and 
other processes that may be affecting the observed line shape. Among the objects that 
have already left the AGB (i.e. pPNs and yPNs) there are only two pA sources. Except 
for these two objects, the rest show evidence of current stellar winds, which supports 
the idea that pAGB winds are generally present in pPNs and yPNs and are most likely 
responsible for the nebular shaping. 

We interpret the peculiar P Cygni like p rofiles observed in some of our o bjects in 



a similar manner as we did for He3-1475 (jSanchez Contreras fc Sahaill200fl ). In this 



scenario the emission and absorption form in two distinct nebular components. The 
broad line emission, characterized by a given intrinsic profile, arises from a compact 
central source. This line emission (and stellar continuum) is scattered by dust in the 
walls of the nebular lobes. The blue-shifted absorption is due to neutral or partially 
ionized outflowing gas inside the lobes absorbing the scattered photons along the line 
of sight. This results in a P Cygni like profile where the blue-shifted absorption is 
produced against the intrinsic emission profile. 

For pE and pcyg sources the FWHM of the intense Ha core emission indicates gas 
motions with velocities in the range [50:200]km s _1 . The mean velocity of the bulk of 
the material producing the blue-shifted absorption in pcyg targets varies between V~50 
and 500km s _1 , however, larger outflow terminal velocities of up to ~ 1000km s _1 
are observed. It is possible that a radial and or latitudinal velocity gradient exists in 
the pAGB outflows and/or that distinct wind components are present. 

Broad Ha emission wings, with widths of up to ±2000km s -1 , are observed in most 
pE and pcyg sources. The yet unclear origin of the very extended wings is inve s tigate d 



following a similar analysis to that performed by lArrieta fc Torres- Peimbertl (120031 ). 
Unlike these authors, we cannot conclude that in our sample Raman scattering is 
the main mechanism for line broadening. The presence of fast stellar winds, which 
is confirmed by P Cygni profiles in more than one third of our targets, is another 
plausible mechanism that could contribute to the broad Ha wings observed. These 



-38- 



winds, however, would have to reach velocities larger than those derived from the P 
Cygni absorption features (i^) to explain the broad Ha wings. 

- The Ha emission filling in the absorption profile of efA indicates moderate speed 
(~50km s _1 ) motions in the nuclear Ha-emitting region. For most efA sources, the 
emission "hump" is displaced bluewards from the narrow stellar absorption feature 
which can be explained by occultation of the receding part of an outflowing stellar 
wind by the central star. 

- We have found differences in the Ha profiles of some of our targets with respect to 
earlier observations, namely, IRAS 19114+0002, IRAS 19475+3119, IRAS 20462+3416, 
and the AGB stars IRC+10216, CIT6, and IKTau. These variations most likely 
represent real changes on the physical properties of the nuclear Ha-emitting region 
(e.g. density, excitation, ionization fraction, geometry and size) presumably induced 
by the evolution of the central star and/or its current stellar wind. 

- We briefly discuss other Hydrogen and non-Hydrogen lines observed towards our tar- 
gets. We note several prominent C 2 Swan bands in the AGB stars IRC+10216 and 
CIT 6 and the pPN IRAS 04296+3429. Some of these bands have been reported for 
the first time in this work. It is also worth mentioning the detection of emission 
by the [Ca II] A7291,7324A doublet in the absorption-line dominated spectrum of 
the YHG IRAS 19114+0002. Observation of these lines, also present in the spectra 
of IRAS 17516-2525, M 1-92, Hen 3-1475, IRAS 22036+5306, IRAS 08005-2356, and 
(tentatively) IRAS 19520+2759, is consistent with dust grain destruction by moderate- 
velocity shocks in the stellar wind. 

- We have estimated the spectral type of the central stars of the objects in our sample 
by comparing their normalized spectra with those of template stars from published 
stellar libraries to investigate the correlation of this fundamental stellar parameter 
with the type of Ha profile. In a number of objects with F- and G-type central stars, 
we found discrepancies (typically of one type) with respect to previous spectral type 
assignations obtained from low spectral resolution studies. Such differences may reflect 
the limitations of low resolution spectroscopy for accurate spectral typing, although a 
real time evolution of the stellar effective temperature cannot be ruled out. 

- The stellar luminosity has been estimated from the luminosity-denpendent O 1 7773A 
infrared triplet, which is observed in absorption in the spectrum of some of our targets 
(13 out of 29). The obtained values, L bo i~5x 1O 3 -1O 5 L , are consistent with pAGB 
stars with initial masses in the range 1-8M , with a significant fraction (10 out of 
13) of objects with masses >3 M . Such massive pAGB objects represent ~30% of 
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our whole sample. The rest of the sources, which is also the majority (~70%), may 
well have low-mass (<3M Q ) progenitors - this may be part of the reason why the O I 
triplet is not observed in absorption. For the highest masses/luminosities derived in 
our sample, the theory of pAGB stellar evolution predicts rapid changes of the effective 
stellar temperature that should be observable in a time-scale of palO-100 years. 

- We investigated correlations between the type of Ha profile and some stellar and en- 
velope parameters. The shape of the Ha line emission is correlated with the stellar 
spectral type as well as the the NIR and IRAS colors. No correlation has been found 
with the chemistry, galactic latitude, or nebular morphology. 

- All sources in which Ha is seen mainly in absorption (i.e., pA and efA) have F-G type 
central stars, whereas sources with intense Ha emission (i.e., pE and pcyg) span a 
larger range of spectral types from O to G, with a relative maximum around B, and 
also including very late C types. The measured equivalent widths of the Ha emission 
in objects with O- and B-type stars are consistent with UV stellar radiation being 
the main ionizing agent. Ha emitters with cooler central stars lack enough ionizing 
radiation, therefore, the emission component in these cases is likely to be formed in the 
de-excitation region behind a shock wave, presumably produced by wind interaction. 

- Pcyg and pE sources are found to exhibit a larger J — K color excess than pA and 
efA objects. Moreover, while the NIR colors of efA+pA sources are consistent with 
the major contribution to the emission in the NIR bands being the reddened stellar 
photosphere, the position of pcyg+pE sources in the (J — H)-(H — K) color-color 
diagram highlights the presence of warm dust. This component of warm dust probably 
results from substantial present-day mass-loss evidenced by the strong Ha emission of 
pE and pcyg targets. 

- Intense Ha emitters (i.e. pcyg+pE) and objects with Ha mainly in absorption (i.e. 
pA+efA) also segregate in the IRAS color-color diagram in a way that the former have 
dust grains with a larger range of temperatures. Such a different temperature range 
may result from: 1) a different distribution of the grain density in the AGB envelope; 
2) a different spatial distribution of grain sizes; 3) a different grain composition; or any 
combination of these situations. 

- The differences found between pE+pcyg and pA+efA sources described above indicate 
dissimilarities in their pAGB and AGB mass-loss histories. The intense Ha emission 
and NIR color excess of pE and pcyg sources may indicate larger pAGB mass-loss rates 
compared with those in efA and pA targets. The lack of a mass/luminosity segregation 
of the two profile sources suggests that pAGB mass-loss is not very dependent on 
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intrinsic stellar properties but may be dictated by extrinsic factors like, for example, 
the presence of a binary companion. The faster winds and earlier spectral types of the 
central stars of pcyg and pE sources suggest that these are older and/or have a larger 
speed traversing the pAGB evolutionary path than pA and efA. 

- We have not found any indications of the extended optical nebulosities (observed in 
most objects) being older for pcyg+pE sources than for efA+pA targets. Such nebu- 
losities show many signs that the jet-sculpting process of the AGB CSE is currently 
active or has been active in the past in all cases. For the only two pA targets in our 
sample, IRAS 19477+2401 and IRAS 17150-3224, the winds that shaped their AGB 
CSEs to its current aspherical morphology are not seen at present. The on-going pAGB 
winds probed by the nuclear Ha emission in pE+pcyg and efA objects, which repre- 
sent the vast majority of our sample, may not be the same that carved and accelerated 
the much more extended (probably older) optical lobes, however, shaping of the inner- 
most layers of the AGB envelope (by interaction with the present-day pAGB winds) is 
probably still at work. 

- In principle, pAGB winds may have been ejected in a continuous or episodic manner 
since the shaping of the AGB envelope began. The fact that pcyg and efA objects with 
central stars of similar late types (F and G) exhibit quite different Ha emission profiles 
may suggest episodic pAGB ejections: pcyg would represent targets that are caught 
at times of intensive jet activity leading to energetic shocks and intense Ha emission, 
whereas efA could be currently in a status of "mild" wind interaction (resulted from 
smaller velocities and/or densities in the pAGB wind?). In pA targets, the pAGB 
winds that presumably shape their aspherical AGB envelopes may have temporarily 
or permanently ended. 
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A. Discrepant stellar spectral type classification 

The stellar spectral type estimated from our ESI data is different from previous assig- 
nations for several objects (§ 14.21) . In the majority of the cases, we find earlier spectral types 
(i.e. larger T eS ) than those reported in the past. The most extreme case is IRAS 19306+1407, 
which has been classified as G5, but displays clear signatures of a B-type star in our ESI 
data (see below). There is also a total of four pPNs with central stars formerly classified as 
G but with ESI spectra most consistent with F types. In two cases the spectral class de- 
rived from this work is later than in previous measurements, namely, IRAS 17317—2743 and 
IRAS 19477+2410, which have gone from a past mid-F to a current G4 and GO classification, 
respectively. Incidentally, we note that these two objects are efA+pA, whereas sources for 
which our spectral type is earlier than in former works are all pcyg+pE. 

Discrepancies between spectral types obtained from low- and high- resolut i on st udies 



has been noticed for oth er pAGB objects and is discussed by, e.g., iDecin et al.l (119981 ) and 



Reddy fc Hrivnaki (119991 ) . Among the objects with discrepant past and our current spec- 
tral type classifications, IRAS 04296+3429 is the only one with additional high-resolution 
spectroscopic line studies in the literature. The fact that the spectral type and T c s for 
IRAS 04296+3429 determined by such studies (including this work) are all in very good 
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agreement but systematically differ from the spectral classification from the low resolution 
spectrum indicates that the changes observed do not correspond to a real time evolution 
of the stellar T e g, but rather reflect the inadequacy/limitations of low-resolution spectra for 
accurate spectral-type determinations. (Note that the previous high-resolution data of this 
object were taken at different epochs, from 1988 to 2003, and the low-resolution spectrum 
was obtained in 1991.) This could also well be the case for objects with low and high- 
resolution spectral types in disagreement by one type, e.g. objects previously classified as F 
(G) versus our current classification as G (F). In other words, we assume that the typical 
error of the spectral type derived from low-resolution spectra can be of one type but unlikely 
larger than that. (Note, however, that a real time variation in T e g cannot be completely 
ruled out.) 

According to this, the spectral type change observed in IRAS 19306+1407 can hardly be 
accounted for by the different spectral resolution. Although it is not impossible, we believe 
that the different spectral type does not result from the rapid evolution of the central star 
towards higher te mperatures. There a re remarkable differences between our spectrum and 



that observed by ISuarez et al.l (120061 ) that suggest that the spectrum reported by these 



authors probably did not correspond to IRAS 19306+1407. This is mainly supported by 



the ab sence of the deep Na I D doublet lines at ~5900A in the spectrum by ISuarez et al. 



(120061 ). The radial velocity derived from the sodium lines in our data, VLSR~30km s , is 
very different from that measured in stellar absorption lines, such as C II6578,6583A and 
He I6678A, which yield ^ LSR ~90-100km s~\ in agreement with the systemic velocity of the 
source derived from CO measurements (Table[3]). The velocity of the sodium lines indicates 
that their origin is most likely interstellar and, therefore, their absence in Suarez et al.'s 
spectrum, which is not explained by their lower S/N, suggests that the object observed by 
these authors is not IRAS 19306+1407. (We note that if the Na ID lines in our spectrum 
were stellar then they should have been even deeper in a G5 star.) Another indication of the 
possible object misidentification is the absence in Suarez et al.'s data of the broad diffuse 
bands at 4430, 5780, 5797, 6284, 6614, 6993 and 7224A observed by us, which may also be 
interstellar. 

In spite of the above reasoning, we would like to note that it is not totally unrealistic to 
expect spectral changes due to rapid pAGB evolution of the central star of IRAS 19306+1407 
in a time scale of ~10 yr. In fact, the high luminosity derived from the deep O I infrared 
triplet indicates a final mass of the central star in the range [>0.84-0.94M Q ], which for 
the standard/empirical initial-final mass relationship, implies an initial mass of Mzamk~5- 



8Mm. According to the theoretical evolutionary models for massive pAGB stars by iBloecker 



(119951 ). a 5-8M star reaches T e ff~20,000K in only 80-10 years, respectively. (For compari- 



son, traversing the same spectral interval would take ~10 3 yr for a Mzams=2M g 
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To the best of our knowledge, evolution towards higher temperatures has been previously 
report ed for two pAGB stars: SAO 85766 (IRAS 18062+2410) and SAO 244567 (Hen 3- 
1357) (jParthasarathy et al.ll2000l . and references therein). For SAO 85766, an increase of the 
temperature from ~8500 K to ~22000 K in less than 25 yr is found, and SAO 244567 has 
turned into a young PN within the last 20-30 years and its central star appears to be evolving 
rapidly into a white dwarf. A very fast evolutio n of the temperature of the central star of 
the YHG IRC+ 10420 has been also reported by iKlochkova et al.l (119971 ) who measured an 
increase of T e g- from ~5700 to ~8600K in only 4 years during 1992 to 1996. 
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Table 1: Journal of observations 



IRAS 


Other 


R.A. (2000) 


Dec. (2000) 


Date observed 


Exposure 


P.A. 


name 


name 


(h m s) 


(°'") 


(yy-mm- 


dd) 


time (s) 


(°) 


01037+1219 


IRC+10011 


01:06:25.98 


+12:35:53.1 


2004- 


-11- 


-08 


1800 


135 


02316+6455 


V656 Cas 


02:35:44.67 


+65:08:58.7 


2004- 


-11- 


-08 


700 


166 


03507+1115 


IK Tau 


03:53:28.87 


+ 11:24:21.7 


2004- 


-11- 


-08 


1250 


120 


04296+3429 




04:32:56.97 


+34:36:12.4 


2004- 


-11- 


-08 


1500 


95 


05506+2414-Sa 




05:53:43.56 


+24:14:44.7 


2004- 


-11- 


-08 


1800 


136 


05506+2414-56 




05:53:43.18 


+24:14:44.1 


2004- 


-11- 


-08 


1200 


45 


08005-2356 




08:02:40.71 


-24:04:42.7 


2004- 


-11- 


-08 


1320 


132 


09452+1330 


IRC+10216, CW Leo 


09:47:57.41 


+13:16:43.6 


2004- 


-11- 


-08 


3600 


22 


10131+3049 


CIT-6, RW LMi 


10:16:02.29 


+30:34:19.1 


2004- 


-11- 


-08 


1470 


20 


17150-3224 


AFGL 6815 


17:18:19.86 


-32:27:21.6 


2003- 


-04- 


-27 


1200 


94+ 


17317-2743 




17:34:53.29 


-27:45:11.5 


2003- 


-06- 


-03 


900 


7 


17423-1755 


Hen 3-1475 


17:45:14.19 


-17:56:46.9 


2003- 


-06- 


-04 


600 


135 


17440-3310 




17:47:22.72 


-33:11:09.3 


2003- 


-04- 


-28 


3600 


115++ 


17441-2411 


Silkworm nebula 


17:47:13.49 


-24:12:51.4 


2003- 


-06- 


-03 


2700 


15 


17516-2525 




17:54:43.35 


-25:26:28.0 


2003- 


-06 


-02 


1920 


164 


18167-1209 




18:19:35.50 


-12:08:08.2 


2003- 


-06- 


-04 


800 


3 


19024+0044 




19:05:02.06 


+00:48:50.9 


2003- 


-06- 


-02 


2400 


138 


19114+0002 


AFGL 2343 


19:13:58.61 


+00:07:31.9 


2003- 


-06- 


-02 


26 


55 


19292+1806 




19:31:25.37 


+ 18:13:10.3 


2004- 


-11- 


-08 


600 


137.4 


19306+1407 




19:32:55.09 


+ 14:13:36.9 


2003- 


-06- 


-04 


1800 


11 


19343+2926 


Ml-92 


19:36:18.91 


+29:32:50.0 


2003- 


-06- 


-02 


1800 


131 


19374+2359 




19:39:35.55 


+24:06:27.1 


2003- 


-06- 


-04 


1800 


15 


19475+3119 




19:49:29.56 


+31:27:16.3 


2003- 


-06- 


-02 


40 


112 


19477+2401 


Cloverleaf nebula 


19:49:54.91 


+24:08:53.3 


2003- 


-06- 


-03 


1800 


91 


19520+2759 




19:54:05.87 


+28:07:40.6 


2003- 


-06- 


-04 


1500 


139 


20028+3910 




20:04:35.98 


+39:18:44.5 


2004- 


-11- 


-08 


600 


165 


20136+1309 




20:16:00.51 


+13:18:56.3 


2003- 


-06- 


-02 


1500 


20 


20462+3416 


LS II +34 26 


20:48:16.64 


+34:27:24.3 


2003- 


-06- 


-02 


1660 


63 


21282+5050 




21:29:58.48 


+51:04:00.3 


2003- 


-06- 


-04 


600 


165 


22036+5306 




22:05:30.28 


+53:21:32.8 


2003- 


-06- 


-04 


2100 


66 


22574+6609 




22:59:18.36 


+66:25:48.3 


2004- 


-11- 


-08 


3900 


130 



tPA changed from 85.9 to 102.1° 
++PA changed from 110.3 to 121.5° 

Note. — All sources were observed with Keck-II/ESI except for IRAS 17150-3224 and IRAS 17440-3310, 
which were observed with Magellan-I/MIKE (see §[2]); Coordinates are 2MASS J2000. 
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Table 2: Stellar and envelope parameters for the different types of Ha profiles 



IH.AS 


Object 


Spectral 


fl2/f25 


f60 


Morphology 1 


c/o' A 


References 


name 


Class 


Type 




(Jy) 


(optical/NIR) 






Pure emission profiles 


(PE) 












09452+1330 


AGB 


C9.5 


2.06 


5652.0 


B vv 


c 


24 


10131+3049* 


AGB 


C4,3e 


2.72 


273.6 


B* 


c 


7, 25 


19374+2359 


yPN 


B3-6 I 


0.24 


70.9 


B 


o 


15, 24, 23 


P Cygni profiles (pcyg) 












05506+2414-Sa 


YSOT 


G9-K2 


0.23 


103.4 


I 


o 


15 16 


08005-2356 


pPN 


F5Ie 


0.35 


29.8 


B* 


c/o 


23, 24 


17423-1755 


yPN 


Be 


0.25 


63.7 




o 


15, 23 


17516-2525 


pPN 


O-B? 


0.45 


100.1 


s 


o 


17, 26 


19024+0044 


pPN 


G0-5 


0.06 


42.5 


M vv 


o 


14, 15 


19306+1407 


yPN 


BO-1 I 


0.06 


31.8 


B* 


c/o 


8, 15 


19343+2926 


pPN 


B2 (+F5) 


0.29 


118.0 


Bw* 


o 


1, 15 


19520+2759 


yPN 


09 I? 


0.39 


207.0 


s 


o 


4, 17 


20462+3416 


yPN 


Bllae 


0.02 


12.1 


E* 


o 


4, 22, 24 


21282+5050 


yPN 


09.5,WC11 


0.68 


33.4 


M * 


c 


2, 15 


22036+5306 


pPN 


F4-7 


0.18 


107.2 


Bw* 


o 


13, 15 


22574+6609 


pPN 


Al-6 I 


0.31 


20.6 


I/Bw 


c 


18 21 


Emission filled 


absorption profiles (efAj 












04296+3429** 


pPN 


F3 I 


0.28 


15.4 


Bw(*?) 


c 


5, 11, 15, 23, 24 


17317-2743 


pPN 


G4I 


0.04 


29.5 


Bw+t 


o 


10, 20, 22 


17440-3310 


pPN 


F3 I 


0.20 


30.1 


B*(w?) 


o 


15, 20 


17441-2411** 


pPN 


F4-5I 


0.22 


106.0 


Bw* 


c? 


11, 15, 22, 23, 24 


18167-1209** 


pPN 


F7I 


0.16 


21.3 


S 


o 


10, 17 


19114+0002 


pPN/YHG? 


F5-7Ia 


0.05 


516.0 


E* 


o 


12, 22, 24, 27 


19475+3119 


pPN 


F3Ib 


0.01 


55.8 


M* 


o 


15, 23 


20028+3910 


pPN 


F3-7 I 


0.20 


143.0 


Bw 


c 


15, 23, 24 


20136+1309 


pPN 


F3-TI 


0.44 


2.1 


S 


O? 


21, 23, 24 


Pure absorption profiles 


(pA) 












17150-3224 


pPN 


F3-7I 


0.18 


268.3 


Bw* 


O 


3, 15, 24 


19477+2401 


pPN 


GO I 


0.20 


27.1 


M 


c 


15, 21, 22 


Ha 


non- detections 














01037+1219 


AGB 


M7 


1.19 


215.2 


E 


O 


7, 9, 15 


02316+6455 


AGB 


M8 


1.53 


45.8 


S 


o 


7, 17 


03507+1115 


AGB 


M6e-M10e 


1.95 


332.1 


I 


O 


7, 9 


19292+1806 


pPN 


B7 


0.10 


28.8: 


Bw 


O 


15, 18 



*Main morphology descriptors defined by Sahai et al. ( 2007bl ). B = bipolar, M = multipolar, E = elongated, 
I = irregular, S = stellar (i.e. unresolved in HST images) , w = central obscuring waist, 7k- = central star evident. 
2 = Oxygen rich; C = Carbon rich; C/O = mixed chemistry. 
"Uncertain profile due to weak Ha emission blended with other lines. 
"""Uncertain profile: could also be pA. 
^YSO serendipitously discovered in our survey of pPNs - see §[T] 

T^No optical image available; our long-slit spectrum shows two nebulosities separated by a dark lane, i.e., most 
likely a bipolar nebula. 



Note. — Boldfaced spectral types have been estimated in this work (see § 14. 2p . 
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Table 3: Parameters of the Ha profile for pcyg and pE sources 









Emission 






Absorption 




IRAS 


Kys(LSR) 


K(LSR) 


FWHM 


FWZI 


w x 


V 


Knax 


FWHM 


name 


(km s^ 1 ) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


(A) 


(km s _1 ) 


(km s _1 ) 


(km s _1 ) 


pcyg srcs 


















08005-2356 


+47 (4) 


+ 101 


190 


-2400 


14 


-85 


-170 


80 


17423-1755 


+48« 


+89 


170 


-4000 


110 


t 


t 


t 


17516-2525 


_1 5 (6) 


+35 


100 


-1400 


105 


-40 


-110 


70 


19024+0044 


+50( 2 ) 


+ 101 


90 


-2400 


40 


-100 


-170 


80 


19306+1407 


+90® 


+ 116 


60 


-2600 


15 


-115 


-170 


50 


19343+2926 


_l(i) 


+19 


170 


-3000 


180 


-530 


-750 


220 


19520+2759 


_1 5 (5) 


-5 


170 


-2600 


130 


-260 


-420 


120 










5) 


5) 


-530 


-800 


250 


20462+3416 


_ 75 (8) 


-43 


80 


-2800 


15 


t 


t 


t 


21282+5050 


+18W 


+17 


50 


-500 


23 


t 


t 


t 


22036+5306 


_ 43 (3) 


+0 


220 


-2500 


30 


-150 


-250 


110 


22574+6609 


_ 64 (D 


-41 


70 


-400 


3 


-50 


-110 


50 


pE srcs 


















19374+2359 


-35® 


-7 


90 


[600:2000] 


4 








09452+1330 


-25® 


-32 


65 


-150 


4 









t Absorption partially masked by emission within the PSF in our ESI spectra (see § 13. lj) 



Note. — References for V sys (LSR): (1) Bujarrabal et al. 2001; (2) Sahai et al. 2005; (3) Sahai et al. 2006; 
(4) Slijkhuis et al. 1991 ; (5) CO data from Sanc hez Contreras et al., in prep.; (6) Nyman et al. 1998; (7) 
Teyssier et al. 2006; (8) llurner k. Drilling! |l984h 
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Table 4: Gaussian analysis of the Ha profile for efA and pA sources 





Narrow absorption 




Broad absorption 




H T > 1 1CG1AT1 
JJjlilloalUIJ 




IRAS 


w x 


FWHM 


y ' A 


AA 


FWHM 


''A 


AA 


FWHM 


name 


(A) 


(km s" 1 ) 


(A) 


(km s _1 ) 


(km s _1 ) 


(A) 


(km s _1 ) 


(km s _1 ) 


17317-2743 


2.22 


140 


2.25 


+ 1(21) 


880 


0.73 


-28(6) 


60 


17440-3310* 


3.94 


70 


0.98 


-10(3) 


200 


4.89 


+4(1) 


100 


19114+0002* 


1.96 


80 








1.56 


+8(3) 


110 


19475+3119 


0.77 


60 


2.22 


-1(3) 


720 


0.60 


-20(2) 


50 


20028+3910 


0.67 


50 


1.27 


-50(17) 


490 


0.32 


-40(10) 


70 


20136+1309 


1.09 


60 


1.10 


-32(3) 


580 


0.79 


-1(2) 


80 


04296+3429 


0.65 


60 


2.04 


+ 10(3) 


640 


0.21 


+7(5) 


125 




0.54 


50 


1.97 


-1(3) 


710 








17441-2411 


0.61 


50 


1.86 


-13(3) 


730 


0.09 


-34(8) 


60 


n 


0.54 


50 


1.88 


-18(3) 


750 








18167-1209 


0.69 


55 


1.94 


-21(2) 


530 


0.16 


-10(3) 


80 




0.56 


50 


1.89 


-27(2) 


520 








17150-3224 


0.50 


50 


1.36 


+30(13) 


530 








19477+2401 


1.09 


110 















*Parameters particularly uncertain (see text § 14. 1|) . 

Note. — Values of AA are relative to the center of the narrow absorption core. AA errors (in km s^ 1 ) are 
indicated in parentheses. 
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Table 5: Equivalent widths of the I 7773A triplet and derived magnitudes 



IRAS 


W x (0 I) 


My 


B.C. 


^bol 


L/D 2 


D 


name 


(A) 


(mag) 


(mag) 


(L Q ) 


(L kpc- 2 ) 


(kpc) 


04296+3429 


1.62 


-5.83 


-0.00 


1.7E+04 


280 


<8 


17441-2411 


1.81 


-6.42 


-0.03 


3.0E+04 


810 


<6 


18167-1209 


1.93 


-6.81 


-0.09 


4.5E+04 


100 


<21 


19024+0044 


1.69 


-6.04 


-0.21 


2.5E+04 


250 


<10 


19114+0002 


3.00 


-10.2 


-0.03 


9.5E+05* 


4950 


<14t 


19306+1407 


0.99 


-3.81 


-2.14 


1.9E+04tt 


320 


<8tt 


19374+2359 


1.98 


-6.95 


-0.95 


l.lE+05 tt 


450 


<13 tf '* 


19475+3119 


2.10 


-7.35 


-0.03 


6.9E+04 


380 


<13 


19477+2401 


1.20 


-4.47 


-0.15 


5.5E+03 


240 


<5* 


20028+3910 


2.04 


-7.15 


-0.03 


5.8E+04 


950 


<8 


20136+1309 


1.26 


-4.67 


-0.03 


5.9E+03 


75 


<9 


20462+3416 


0.67 


-2.81 


-2.14 


7.4E+03 


130 


<8 


22574+6609 


2.10 


-7.34 


-0.13 


7.6E+04tt 


170 


<2ltt 



^Uncertain, outside of the My applicability domain (see § 14. 3|) . 
tt Uncertain, outside of the spectral type applicability domain (see § 14. 3p . 
* Interstellar Ay expected to be very high. 



Fig. Set 1. Observed spect 
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Fig. 1. — The spectrum of IRAS 19306+1407 observed with ESI is shown here as an example. 
In some panels, the spectrum is also plotted using a smaller flux scale for clarity (grey line). 
Spectra for all sources are available in the electronic edition of the journal. 
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Fig. 2. — Spectra in the region around the Ha line for objects with Ha emission, i.e. pcyg 
and pE sources (TableEJ). In some cases, the spectrum is also plotted using a smaller scale 
for clarity (red line). 
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BUT™ 




Fig. 3. — Spectra in the region around the Ha line for objects with Ha absorption, i.e. efA 
and pA sources (Table^ . The spectra of IRAS 17150-3224 and IRAS 19477+2401 have 
been boxcar smoothed (using a flat-topped kernel of 5 pixels) to increase S/N. 
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Fig. 4. — Spectra in the region around H/3 for sources with Ha P Cygni like profiles 
(Tabled . The spectra of the pcyg sources IRAS 17516-2525 and IRAS 22574+6609, for 
which signal-to-noise ratio is <3 in this region are not plotted. For IRAS 20462+3416 and 
IRAS 21282+5050, the spectrum is also plotted using a smaller scale for clarity. 
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Fig. 5. — Fit to the Ha wings of pcyg sources by a function of the type I a oc A~ 2 (dashed line), 
which is expected for Raman scattering (§ 15.11) . The broad line wings of IRAS 19306+1407 
and IRAS 20462+3416 cannot be reproduced by the fit. 
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Fig. 6. — Fit to the Ha-profile for efA sources using Gaussian functions (dotted line= 
emission component; dashed line= absorption components; solid line= resultant profile). 
Line parameters obtained from the fit are given in TableHl For IRAS 19114+0002, one 
gaussian function has been fitted to the Ti I line adjacent to Ha to better reproduce the 
observed profile. 
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Fig. 7. — Gaussian fit to the Ha-profile for pA sources (dashed line= absorption components; 
solid line= resultant profile). Parameters of the fit are given in TableHJ 
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Fig. 8. — Positions of our targets with the infrarred O I triplet seen in absorption (Ta- 
ble|5]), in the HR diagram (filled squares: pcyg+pE; open squares: efA+pA). Errorbars 
of Alog(T ef j)=0.05 and Alog(L/L )=O.5 have been adopted. The luminosities derived for 
objects with B-type central stars and for the YHG IRAS 19114+0002 are particularly un- 
certain since these objects could be outside of the applicability domain of the My-W\(0 I) 
relationship (see text § 14.31) . 




Fig. 9. — SEDs of the objects with infrared I triplet absorption (TableEJ)- Big circles rep- 
resent GSC2, 2MASS, MSX and IRAS data points; Small circles are used for ISO and IRAS 
LSR spectra. Two blackbody curves (dashed lines) have been fitted to the reddened stellar 
photosphere and cool dust emission components. The two-blackbody SED model corrected 
for intestellar extinction is also shown (green solid line); an typical value of Ay =2. 5 mag is 
used for reddening correction. 
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OBAFGKMC 
Spectral Type 

Fig. 10. — Histogram showing the distribution of spectral types for pcyg+pE (grey-filled 
bars) and efA+pA (blank bars) sources. 
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Fig. 11. — {Top) Distribution of our targets in the (J-H)-(H-K) color-color diagram (filled 
squares: pcyg+pE; open squares: efA+pA; crosses: non-Ha detections). The reddening 
vector for v4j=lmag is indicated. The solid line represents the locus of blackbody emitters 
for temperatures in the range 6000-750 K; the tick marks correspond to the temperatures 
indicated. (Bottom) Histograms of the J — K and H — K colors for pcyg+pE and efA+pA 
sources (same color code as in Fig. fTOl . 



-66- 



6 - 



~i — 1 — i — 1 — r 

□ efA+pA 
*non— Ha 



i — | — r 



"i — 1 — i — 1 — i — 1 — r 



4 - 



I 



- 



T 



^1 




■ 



C— stars: 



lass III 
class I 



J i L 



00 BO AO FO GO KO MO CO CIO 
Spectral Type 



Fig. 12. — J — K vs. spectral type diagram showing the observed colors for our targets 
(symbols) and th e intrinsic colors for luminosity classes I, III, and V, and for C-rich stars 
(solid lines; from IDucati et all (120011 )). Horizontal errorbars represent uncertainties in the 
spectral classification of the central stars. For the J — K color, an error of O.lmag has been 
adopted. Symbol code as in previous figures. 
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Fig. 13.— (Left) Distribution of our targets in the IRAS [12] -[25] vs. [25] -[60] color 
diagram. The solid line represents the blackbody colors for temperatures in the range 100- 
860 K. The dashed line represents the colors for optically thin dust emission for the same 
range of temperatures and dust optical depth r u oc (see § 14. 4p . (Right) [12] — [25] histogram 
— same color-symbol code as in Fig.fTTl 



- 68 - 



-i 1 1 1 1 1 r 



"i r 



~T 1 — fj 



50 



a> 



■ 

□ 




■ X 



□ 



□ 



X 




Fig. 14. — Distribution of our targets in the Galaxy. Symbol code as in previous figures. 



